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FOREWORD 

This  report  is  an  adaptation  of  Don  A.  Webster's  thesis, 

"Saturation  of  Plane  Acoustic  Waves  and  Notes  on  the  Propagation  of 

Finite-Amplitude  Spherical  Waves, " vhich  was  written  for  the  degree 

of  Master  of  Science  in  Engineering  at  The  University  of  Texas  at 

Austin.  Mr.  Webster  was  enrolled  in  the  Department  of  Electrical 

Engineering,  and  his  degree  was  granted  in  December  1976. 

The  research  was  carried  out  at  the  Applied  Research  Laboratories 

and  was  supported  by  the  Air  Force  Office  of  Scientific  Research  under 

Contract  F44620-76-C-0040,  by  the  National  Aeronautics  and  Space 

Administration  under  Contract  NASl-l4l6o,  and  by  the  National  Oceanic 

uJ 

and  Atmospheric  Administration  under  Grant  04-5-022-12.  Technical 
monitors  were  Lt.  Col.  R.  C.  Smith  and  Lt.  Col.  L.  W.  Ormand  for  AFOSR, 
Dr.  J.  M.  Seiner  for  NASA,  and  Dr.  F.  F.  Hall  for  NOAA. 

Readers  particularly  Interested  in  the  outdoor  propagation 
measurements  reported  briefly  in  Appendices  B and  C may  wish  to  consult 
Mark  A.  Theobald's  thesis,  "Experimental  Study  of  Outdoor  Propagation 
of  Spherically  Spreading  Periodic  Acoustic  Waves  of  Finite  Amplitude." 
This  thesis,  which  is  Ref.  28  in  the  present  work,  is  soon  to  be  issued 
as  Applied  Research  Laboratories  Technical.  Report  ARL-TR-77-5  (January 
1977)  and  subsequently,  in  somewhat  different  form,  as  a NASA  Contractor 
Report.  Experimental  details  as  well  as  several  additional  outdoor 
experiments  are  given  in  Theobald's  thesis. 


David  T.  Blacks tock 
Supervisor 


! ABSTRACT 
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/Tte  saturation  of  plane  acoustic  waves  In  an  air-filled  tube.-ts-e- 
inreabigated.?  Saturation^jiihieh-t?  a state  in  which  the  signal  amplitude 
at  a field  point  approaches  a limiting  value,  independent  of  source  ampli- 

Cjfitv  OlU  fJctrx 

tudej  ~is  caused  by  nonlinear  propagation  distortion.  Wr  observed  it  in  a 
50-m  long,  5-cm  i.d.  progressive  wave  tube.  Measured  saturation  curves 
were  obtained  for  the  following  source  conditions:  sound  pressure  levels 
from  110  to  163  dB,  frequencies  from  500  Hz  to  3500  Hz.  The  model  equa- 
tions of  Merklinger  and  Rudnick  were  adapted  and  solved  to  give  theoretical 
predictions  for  the  fundamental  component  of  the  wave.  The  Merklinger 
model  is  found  to  give  good  results  at  low  levels,  and  the  Rudnick  model 

at  high  levels.  measurements  generally  confirm  fehg^^eoretical 

( efCicf  c-P- 

predictions^  Aspersion  due  to  the  tube  wall  boundary  layer,  is  not 
included  in  either  of  the  theoretical  models.  Its  effect '•was  analyzed 
separately  and  found  to  be  not  great  for  the  fundamental  component. 

Certain  aspects  of  the  propagation  of  spherically  spreading 
waves  of  finite  amplitude  were  also  studied.  A graphical  method  for 
assessing  the  importance  of  nonlinearity  on  spherical  waves  is  discussed. 

New  theoretical  results  for  weak  waves  and  for  strong  waves  were  found. 

In  the  first  case  a high  order  perturbation  solution  of  Burgers'  equation 

«■  — ^ 

was  obtained.  In  the'  second,  the  Rudnick  model  equation  for  the  funda-  

v/'j>  ?-  ’’a  f 

mental  component  of  a sawtooth  spherical  wave  was  solved.  In  both  cases,  i.*t 
measurements  taken  in  two  outdoor  propagation  experiments  confirmed  the  


theoretical  results. 
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CHAPTER  1 


INTRODUCTION 


In  linear  acoustic  theory  the  amplitude  of  the  received  signal 
at  a fixed  point  varies  in  direct  proportion  to  the  source  amplitude.  As 
the  source  amplitude  is  increased,  however,  losses  due  to  finite  amplitude 
effects  develop.  The  received  level  increases  less  rapidly  than  the 

source  level.  Eventually  the  finite  amplitude  losses,  or  extra  attenua- 

# 

tion,  become  so  severe  that  the  received  signal  approaches  an  upper 
limit,  independent  of  source  amplitude.  The  wave  is  said  to  have 
saturated. 

The  terminology  used  to  describe  the  development  of  acoustic 
saturation  may  be  made  clear  by  Fig.  1-1.  What  is  plotted  here  is  the 
received  sound  pressure  level  (SPL)  at  a fixed  distance  versus  the  source 
SPL.  Such  a plot  is  commonly  referred  to  as  an  amplitude  response  or 
input-output  curve.  At  sufficiently  low  source  levels  the  wave  propagates 
as  a small  signal  and  the  input-output  curve  (solid  curve)  has  unit  slope. 
Eventually  the  curve  begins  to  bend  over  and  depart  from  the  extrapolated 
small-signal  curve  (dashed  line),  signifying  that  the  approach  to  satu- 
ration has  begun.  The  extra  attenuation  in  decibels,  conmonly  called 
EXDB,  is  the  difference  between  the  extrapolated  small-signal  curve  and 
the  input-output  curve.  The  saturation  level  is  the  asymptotic  value  of 
the  received  level  in  the  limit  of  high  source  levels.  In  practice,  of 


Extra  attenuation  here  is  defined  as  attenuation  over  and  above  that 
observed  for  a small-signal  wave. 


1 


course,  the  asymptote  can  never  be  attained.  It  is  therefore  useful  to 
define  a point  at  which  the  received  level  differs  from  the  saturation 
level  by  some  prescribed  amount  (for  example,  1 dB).  This  point  defines 
the  beginning  of  the  plateau  region,  where  the  input-output  curve  i6 
nearly  horizontal.  In  thie  region  the  wave  is  practically  saturated. 


FIGURE  1-1 

THE  DEVELOPMENT  OF  SATURATION 

The  extra  attenuation,  which  is  responsible  for  saturation,  is 
caused  by  a transfer  of  the  wave's  energy  to  higher  frequency  components 
where  it  is  more  efficiently  dissipated.  The  energy  transfer  may  be 
traced  to  the  dependence  of  the  propagation  speed  dx/dt  of  the  wave  on  the 
local  particle  velocity  u.  The  propagation  speed  of  a particular  wavelet 
(hence  a particular  value  of  u)  is* 


(1-1) 


u=constant  co  + 

Here  cq  is  the  small-signal  sound  speed  and  p is  the  parameter  of 
nonlinearity  (p=(7+l)/2  for  gase3,  where  7 is  the  ratio  of  specific  heats). 
The  dependence  of  dx/dt  on  u is  a consequence  of  two  phenomena.  First, 
since  sound  propagates  longitudinally,  the  wave  is  carried  along  by  the 
motion  of  the  medium  (i.e.,  the  sound  wave  is  convected).  Second,  because 
of  the  nonlinearity  of  the  pressure-density  relation,  wavelets  associated 
with  the  compression  phase  of  the  wave  travel  faster  than  wavelets  asso- 
ciated with  the  rarefaction  phase.  These  two  effects  reinforce  each 
other  and  cause  a cumulative  distortion  of  the  wave  as  it  propagates. 

For  example,  in  the  absence  of  ordinary  dissipation,  a sinusoidal  wave 
will  distort  into  a sawtooth  wave.  Hie  harmonic  distortion  necessary  to 
produce  the  sawtooth  comes  at  the  expense  of  the  fundamental  component. 

In  other  words,  energy  transfer  takes  place. 

Most  of  the  work  that  has  been  done  on  acoustic  saturation 

2 

concerns  spherically  spreading  sound  waves.  Allen,  who  made  measurements 
on  spherical  waves  in  air,  reported  the  first  experimental  evidence  of 
acoustic  saturation  in  1950*  Allen's  measurements  of  the  first  six 
harmonics  of  the  received  signal  show  a definite  approach  to  saturation. 

At  a distance  of  2 m from  the  source  Allen  found  that  an  Increase  in 
source  SPL  from  158  dB  to  l6l  dB  yielded  only  a 0.5  dB  increase  in  the 
received  level.  An  experiment  on  saturation  of  spherical  waves  in  water 
is  given  in  Ref.  3»  this  reference  also  contains  an  excellent  historical 

review  of  other  work  relating  to  acoustic  saturation  of  spherical  waves. 

* 

In  the  remainder  of  this  work  unless  otherwise  noted  the  SPL  reference 
pressure  is  0.0002  pbar. 


dx 

dt 
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Little  experimental  work  has  been  done  on  saturation  of  plane 
waves.  Plane  waves  are  usually  achieved  either  by  generating  a collimated 
beam  in  an  open  medium  or  by  confining  the  sound  wave  in  a tube.  Lester 
made  measurements  on  a collimated  ultrasonic  beam  in  water.  His  measure- 
ments confirm  the  plateau  region  shown  in  Fig.  1-1  though  the  plateau 
turnB  down  unexpectedly  at  the  highest  amplitudes.^  Although  much  work  on 
plane  finite-amplitude  waves  has  been  done  in  tubes  (see,  for  example. 

Refs.  5-9),  a specific  investigation  of  acoustic  saturation  has  not  been 

9 

made.  Pernet  and  Payne,  for  example,  made  measurements  of  intense  sound 
in  tubes,  but  they,  as  many  others,  were  primarily  interested  in  the 
harmonic  growth  with  distance  and  hence  did  not  investigate  saturation 
explicitly. 

The  object  of  this  research  is  to  investigate  saturation  of 

plane  acoustic  waves  in  tubes.  Several  model  theories  for  the  fundamental 

component  of  the  wave  are  discussed.  Solutions  of  the  differential 

equations  based  on  these  models  are  compared  with  data  taken  in  a plane 

q 

wave  tube  in  air.  A numerical  algorithm  developed  by  Pestorius  to  com- 
pute the  wave  shape  at  any  propagation  distance  is  also  compared  with  data 
from  the  plane  wave  tube. 

Although  we  have  been  mainly  concerned  in  this  research  with 
plane  waves,  some  new  results  have  been  obtained  for  spherical  waves. 

These  results  are  given  in  the  appendices.  In  Appendix  A a graphical 
method  for  assessing  the  importance  of  nonlinear  effects  on  spherical 
waves  is  described.  This  method,  which  is  a simplification  of  one 
described  in  Ref.  3,  allows  one  to  classify  nonlinear  effects  as  weak, 
moderate,  or  strong.  A perturbation  solution  useful  for  weak  waves  is 
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presented  in  Appendix  B.  A solution  valid  for  strong  waves  is  presented 
in  Appendix  C.  Both  the  weak  and  strong  wave  solutions  are  compared  with 
data  from  freefield  propagation  experiments. 


CHAPTER  2 


REVIEW  OF  THE  LITERATURE 


A.  Introduction 

In  this  chapter  we  review  several  of  the  theoretical  methods 
which  may  be  applied  to  the  problem  of  propagation  of  plane  acoustic  waves 
of  finite  amplitude.  Our  goal  is  to  find  an  expression  for  the  funda- 
mental component  of  an  initially  sinusoidal  wave.  Both  nonlinear  effects 
and  ordinary  dissipation  are  taken  into  account.  The  specific  problem  we 
consider  is  the  propagation  of  progressive  waves  contained  in  a gas-filled 
tube.  The  analysis  in  Section  C of  this  chapter  shows  that  for  our 
experimental  conditions  only  boundary  layer  absorption  is  important.  We 
shall,  however,  in  two  cases  present  results  valid  only  for  a thermo- 
viscous  gas.  These  results  are  included  because  the  methods  used  may,  in 
some  cases,  be  recast  in  terms  of  boundary  layer  absorption. 

B.  Review  of  Weak-Shock  Theory 

Because  extensive  use  is  r-ade  here  of  weak-shock  theory,  an 
outline  of  its  basic  features  is  appropriate.  The  exact  one-dimensional 
equation  for  plane  progressive  waves  in  an  ideal  lossless  gas  is1 


Given  the  boundary  condition 

u(0,t)  = g(t)  , (2-2a) 

Poisson10  found  the  solution  (generalized  here  for  an  adiabatic  gas)  to  be 
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An  extension  of  the  Poisson  solution  (Eq.  2-3),  weak-shock  theory  enables 
us  to  calculate  the  wave  shape  when  shocks  are  present.  Through  the  use 
of  the  Rankine-Hugoniot  shock  relations,  which  state  conservation  law6, 
shock  dissipation  is  included.  An  excellent  account  of  weak-shock  theory 
is  given  by  Blackstock.  We  give  here  a brief  account  of  the  solution 
corresponding  to  a sawtooth  wave  shape.  The  inclusion  of  this  solution 
is  essential  here,  since  it  is  the  formation  of  the  sawtooth  wave  shape 
which  ultimately  leads  to  saturation.  We  consider  a source  excitation  of 
the  form 

p(0,t)  = p1Q  sin  cut  . (2-3) 


Here  p,  u>,  and  p^  are  the  acoustic  pressure,  angular  frequency,  and 
pressure  amplitude,  respectively,  of  the  source.  A discontinuity  in  the 
wave  (a  shock)  is  first  formed  at  a distance  x given  by12 


x = 


1 

Pek 


(2-4) 


Here  e=p^^ P0C0^>  PQ  1b  the  static  density,  and  k=cu/cQ  is  the  wave  number. 

An  approximate  solution  of  Eqs.  2-1  and  2-3  valid  in  the  more  remote 
— 12 

region  x>3x  is 


HO  Vs*  i 

P = Y~rz  2*  n Bin  na*'  • ( 2- 5) 

where  a=x/x.  Equation  2-5  is  the  Fourier  series  representation  for  a 
sawtooth  wave  of  amplitude  xp1(^(l+o).  The  distance  x=3x  (a=3)  is 
conmonly  called  the  sawtooth  distance,  and  the  region  of  validity  of 


Eq.  2-5  is  called  the  sawtooth  region. 
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As  will  be  shown  in  Chapter  3>  Section  B,  Eq.  2-5  may  be  used  to 
indicate  the  development  of  saturation.  In  many  cases,  however,  ordinary 
absorption  must  also  be  taken  into  account.  In  the  following  section  we 
examine  the  role  of  ordinary  absorption  and  dispersion  in  our  experiments. 


C.  The  Role  of  Ordinary  Absorption  and  Dispersion 

Ordinary  absorption  may  be  expected  to  be  important  in  our 

0 

measurements.  Indeed,  Pestoriu6  found  it  necessary  to  include  both 
absorption  and  dispersion  in  his  analysis  of  finite-amplitude  propagation 
in  tubes.  Some  of  the  effects  of  absorption  and  dispersion  on  a small- 
signal  wave  are  discussed  here  so  that  the  role  of  each  in  the  nonlinear 


case  may  be  better  understood. 

The  problem  of  absorption  and  dispersion  of  plane  sound  waves  in 

15  * 

tubes  was  first  solved  by  Kirchhoff.  For  a wide  tube  one  finds  that  the 


absorption  coefficient  a and  phase  velocity  c 


ph 


are  given  by 


and 


'ph 


1 + c 


w° 


(2-6) 


(2-7) 


* lk 

A wide  tube  is  a tube  for  which  the  boundary  layer  thickness  A is  small 

compared  to  a.  We  define  the  boundary  layer  thickness  A as  the  distance 
from  the  tube  wall  to  the  point  where  the  particle  velocity  amplitude  has 
reached  a value  equal  to  1/e  times  its  mainstream  value.  It  turns  out 
that  A=  yj  2 7/ 05.  The  requirement  A«a  thus  Imposes  a lower  limit  on  the 
frequencies  for  which  the  wide  tube  formulas  may  be  used.  The  tube  must 
not,  however,  be  so  vide  that  losses  in  the  fluid  mainstream  are  compara- 
ble to  the  boundary  layer  losses.  The  latter  requirement  may  be  stated  as 
a«l/ k^A.  15  Hence  a vide  tube  must  meet  the  requirement  that  A«a«l/k^A. 
For  example,  to  be  classified  wide  over  the  frequency  region  100  Hz  to 
30  kHz,  the  tube  must  have  a radius  in  the  range  0.02  cm  « a « 30  cm. 

The  pipe  we  use  for  our  experiments  (see  Chapter  4)  has  a radius  of  2.5  cm. 
Our  pipe,  therefore,  is  a wide  pipe  from  100  Hz  to  30  kHz. 
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Here  a is  the  tube  radius,  v 1b  the  kinematic  viscosity,  and  Pr  Is  the 

Prandtl  number.  From  Eqs.  2-6  and  2-7  ve  see  that  c . decreases  with 

ph 

frequency,  asymptotically  approaching  co  in  the  high  frequency  limit. 

The  absorption  coefficient  a varies  aB  the  square  root  of  oj.  From  these 
observations  we  conclude  that  a nonsinusoidal  signal  propagating  in  a tube 
will  suffer  both  phase  and  amplitude  distortion. 

The  importance  of  ordinary  absorption  in  our  experiments  is 
easily  seen  from  Eq.  2-6.  The  tube  we  use  for  our  measurements  has  a 
2.5  cm  radius,  and  the  (fundamental)  frequency  range  of  interest  is  500  Hz 
to  4000  Hz.  At  a frequency  of  5.5  kHz,  for  example,  the  absorption 
coefficient  0!  is  about  0.072  Np/m.  Therefore,  for  the  available  measure- 
ment distance,  26  m,  we  could  expect  16.5  dB  of  small-signal  attenuation 
alone.  We  conclude  that  a reasonable  theory  for  our  experiments  must 
include  the  effect  of  ordinary  absorption. 

The  importance  of  dispersion  may  be  seen  from  an  examination  of 
the  phase  distortion  suffered  by  a nonsinusoidal  signal.  Since  we  are 
concerned  with  saturation  in  our  experimental  work,  the  main  waveform  of 
interest  besides  the  sine  wave  is  the  sawtooth.  McKtttrick  et  al.1^  were 
able  to  explain  asymmetric  waveforms  observed  in  a high  intensity  tube 
experiment  by  considering  the  propagation  of  a small-signal  sawtooth  wave. 
We  present  an  analysis  similar  to  theirs  here  because  it  makes  clear  the 
effect  of  absorption  and  dispersion  on  the  overall  wave  shape.  We  assume 
a sawtooth  wave  of  peak-to-peak  amplitude  pQ  exists  at  x=0.  The  nor- 
malized acoustic  pressure  p/pQ  may  then  be  expressed  as  the  Fourier  series 

KL  " « £ » ,1° 


(2-8) 


is  given  at  distance  x by' 


If  Eqs.  2-6  and  2-7  are  applied,  wPQ 
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(2-9) 


Figure  2-2  shows  the  sum  of  the  first  one  hundred  terms  of  Eq.  2-9  for 
three  cases.  The  top  waveform  is  one  in  which  there  is  no  absorption  or 
dispersion,  i.e.,  ax  has  been  set  to  zero.  Note  the  Gibbs'  overshoot  at 
both  peak  and  trough.  For  the  middle  and  bottom  waveforms  the  value  of  ax 
is  0.5.  In  the  middle  waveform  absorption  is  included  but  not  dispersion, 
i.e.,  the  phase  angles  [iixx  are  ignored.  The  result  is  that  the  waveform 
is  rounded  symmetrically  at  both  peak  and  trough.  In  the  bottom  waveform  " 
both  absorption  and  dispersion  are  included.  The  wave  asymmetry  is  clear; 
the  peak  is  rounded  while  the  trough  remains  sharp. 

We  now  ask  what  effect  dispersion  will  have  on  a wave  of 
finite-amplitude.  Dispersion  clearly  alters  the  wave  shape  in  the  linear 
case  and,  as  McKittrick  et  al.  found,  causes  the  waveforms  to  be  asymmet- 
ric in  the  finite -amplitude  case  as  well.  What  we  are  mainly  concerned 
with,  however,  are  the  harmonic  amplitudes.  For  a small-signal  wave, 
dispersion  does  not  alter  the  harmonic  amplitudes,  i.e.,  the  spectra  of 
the  middle  and  bottom  waveforms  in  Fig.  2-1  are  identical.  The  propaga- 
tion of  a finite-amplitude  wave,  however,  is  highly  dependent  upon  the 
8 

wave  shape.  Consider,  for  example,  the  propagation  of  two  waves,  one  a 
sawtooth  and  the  other  an  inverted  sawtooth.  Although  both  have  the  same 
initial  spectra,  the  finite-amplitude  distortion  of  the  two  waves  is 
vastly  different.  The  reason  is  that  the  phases  of  the  Fourier  components, 
and  therefore  the  wave  shapes,  are  different.  This  is  admittedly  a rather 
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SAWTOOTH  WAVE 


ABSORPTION  ONLY 


BOTH  ABSORPTION 
AND  DISPERSION 


FIGURE  2-2 

EFFECTS  OF  ABSORPTION  AND  DISPERSION 
ON  A SMALL-SIGNAL  WAVE 


, 
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drastic  example,  given  merely  to  point  out  that  dispersion  will  probably 
have  some  effect  on  the  propagation  of  the  wave.  A qualitative  discussion 
of  the  effect  of  dispersion  is  given  in  Chapter  6.  The  discussion  given 
there  is  based  on  the  bottom  waveform  of  Fig.  2-2. 

The  effects  of  both  dispersion  and  ordinary  absorption  are 
usually  included  in  a theoretical  analysis  in  one  of  two  ways,  both  of 
which  are  treated  briefly  in  Section  D.  In  the  first  method  the  effects 
are  included  as  specific  terms  in  the  differential  or  integral  equation 
governing  the  wave  motion.  Unfortunately,  however,  the  equation  is 
generally  very  difficult  to  solve.  Another  method  is  to  include  the 
effects  stepwise  in  a computer  algorithm  that  is  based  on  a mathematical 
model  of  the  wave  distortion  process.  Such  an  algorithm  has  been  described 

Q 

by  Pestorius.  Although  the  algorithm  is  not  efficient  for  computing 
saturation  curves,  we  do  make  use  of  it  in  Chapter  6.  There  a quantita- 
tive discussion  of  the  effects  of  dispersion  are  given. 

A completely  different  approach  is  to  ignore  dispersion  entirely. 
Several  ad  hoc  models  have  been  proposed  in  which  ordinary  absorption  is 
accounted  for  but  not  dispersion.  Perhaps  surprisingly,  several  of  these 
models  have  enjoyed  considerable  success.  Appropriate  ones  Eire  discussed 
in  Section  D. 

D.  Contributions  of  Other  Researchers 

1.  A Computational  Method  - Pestorius1  Algorithm 

In  a computer  algorithm  developed  by  Pestorius,  explicit  account 

Q 

is  taken  of  nonlinear  effects,  ordinEury  absorption,  and  dispersion.  In 
this  algorithm  the  source  wave  is  "computer  propagated"  a short  distance 


Ik 


Ax;  the  basis  of  the  computation  is  weak-shock  theory.  The  resulting  time 
waveform  is  then  transformed  to  the  frequency  domain,  where  corrections 
are  made  for  small-signal  absorption  and  dispersion.  After  a transfer 
back  to  the  time  domain,  another  short  propagation  step  is  made.  In  this 
way  the  time  waveform  may  be  computed  for  any  propagation  distance.  The 
validity  of  this  algorithm  has  been  verified  by  Pestorius  for  both 

Q 

initially  sinusoidal  waves  and  noise.  As  was  mentioned  in  Section  C,  we 
make  use  of  this  algorithm  in  Chapter  6 to  find  the  effect  of  dispersion 
on  the  harmonic  amplitudes.  The  algorithm  was  not,  however,  used  to 
compute  amplitude  response  curves  because  a separate  computer  run  would 
have  been  required  for  each  increment  in  source  SPL. 


2.  Burgers'  Equation 

a.  Mainstream  Absorption 
17 

Lighthill  reduced  the  equations  of  motion  for  progressive 

plane  waves  in  an  unbounded  thermoviscous  gas  by  a series  of  approxima- 

18 

tions  to  Burgers'  equation.  Mendousse  cast  Burgers'  equation  in  the 
following  form,  which  is  equally  valid  but  better  Buited  for  boundary 
value  problems: 


5 c2 


o 


u 


du 

3F 


a 


0} 


(2-10) 


where  a is  the  thermoviscous  absorption  coefficient  at  angular  frequency 
15 

a).  It  will  be  seen  that  the  left-hand  side  of  Eq.  2-10  is  an  approxi- 
mate form  of  Eq.  2-1.  The  right-hand  side  is  the  dissipation  term 
appropriate  for  mainstream  absorption.  Though  exact  solutions  of  Eq.  2-10 
are  known,  they  are  not  relevant  for  our  problem  because  in  our  case  the 
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attenuation  is  due  to  boundary  layer  effects,  not  mainstream  effects.  The 
spherical  wave  version  of  Eq.  2-10  is,  however,  treated  in  Appendix  B by 
a perturbation  method. 


b.  Boundary  Layer  Absorption 

Blackstock1  found  a "Burgers -like"  equation  for  finite-amplitude 
waves  subject  to  boundary  layer  attenuation. 


iE  _ _g_  „ , fT  /*“  d u(x.t'-\)  d\ 

^ c2  " St ' 7= 


(2-11) 


Here  a is  the  boundary  layer  absorption  coefficient  at  angular  frequency 
see  Eq.  2-6.  Solutions  of  Eq.  2-11  would  be  ideal  for  our  work. 
Unfortunately,  however,  only  low-order  perturbation  solutions  have  been 
found  and  they  are  not  applicable  in  our  study. Evaluation  of  Eq.  2-11 
by  numerical  means  was  not  pursued.  Coppens  and  Sanders20  also  found  a 
"Burgers-like"  equation  for  tubes.  Their  equation  is,  however,  less 
general  than  Eq.  2-11  in  that  it  is  limited  to  periodic  waves. 


3*  Ad  Hoc  Models 

A great  body  of  the  literature  regarding  plane  propagation  of 
finite-amplitude  waves  in  a lossy  medium  is  devoted  to  the  subject  of 
ad  hoc  models.  The  models  presented  here  illustrate  some  of  the  methods 
of  blending  nonlinear  effects  and  ordinary  dissipation.  No  one  has  yet, 
unfortunately,  found  a way  to  include  the  effects  of  dispersion. 
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a.  Merklinger 's  Method 
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Merklinger  applied  Westervelt 1 s equation 


— = - 20^1 l - (p(x,t' )q(x,t')) 


to  the  problem  of  plane  wave  propagation  in  a viscous  fluid.  Here  is 
the  intensity  of  the  fundamental  frequency  component  at  distance  x,  a».d 
is  the  absorption  coefficient  (appropriate  for  a viscous  fluid)  at  the 
fundamental  frequency.  The  angled  brackets  indicate  a time  average  over 
one  period  of  the  wave,  and  q is  the  source  function  given  by 


= —§■  n 4*  Ipfat)  ‘ 

r\  ‘ 


P c 
o o 


Merklinger  assumed  that  the  nonlinear  decay  of  the  fundamental  pressure 
component  p1(x,t')  is  primarily  due  to  work  that  it  does  on  the  second 
harmonic  pressure  component  p2(x,t')-  That  is,  Merklinger  approximated 


p = p^Xjt')  + p2(x,t')  . (2-1* 

For  p2(x,t'),  Merklinger  used  a second  order  perturbation  solution^  of 
the  nonlinear  wave  equation  for  a viscous  fluid.  The  boundary  condition 
is  given  by  Eq.  2-3. 


2q  -2a, x 

, ..X  pi0^  e 1 

P2(x,t')  = T 

O O 


-to^x 


sin  2cut' 


We  here  define  Pjj  to  be  the  pressure  amplitude  of  the  nth  harmonic.  If 
the  time  dependence  is  to  be  included,  we  shall  use  the  form  pn(x,t'). 


r 


Writing  p^Xjt')  as 


P1(x,t')  = p-^x)  sin  cut’ 


(2-16) 


2 ~2a±x 


and  reinterpreting  P10  e in  Eq.  2-15  as  P-j^  (x),  we  may  write  p2(x,t') 


P2(x,t‘)  = px2(x) 


4p  C 

o o 


s;  (■  • •*“) 


sin  2<iit ' . (2-17) 


The  differential  equation  for  p1  may  be  found  from  Eqs.  2-12  through  2-17. 


The  result  is 
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(2-18) 


The  solution  of  Eq.  2-18  consistent  with  the  boundary  condition  p^(0)=p^ 


P10  e 


jZ  / -ai  x 


(2-19) 


where  r=l/a^x.  The  assumptions  inheieut  in  Eqs.  2-14  and  2-15  should  be 
expected  to  limit  the  validity  of  Eq.  2-19.  13118  problem  is  addressed  in 
Chapter  Section  B.  We  note  in  passing  that  Merkllnger  has  found 

Eq.  2-19  bo  be  in  good  agreement  with  an  exact  solution  of  Burgers' 
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equation  for  a wide  range  of  source  amplitudes. 

In  using  Eq.  2-15  for  p2(x,t')  we  have  assumed  a viscous  fluid. 

To  make  use  of  Merkllnger ’s  method  for  the  case  of  boundary  layer  absorp- 
tion we  must  find  an  appropriate  expression  for  p2(x,t' ).  Ihuras,  Jenkins, 


18 


and  O'Neil^  derived  an  equation  in  which  no  assumption  is  made  of  a 
specific  absorption  mechanism.  Their  method  is  presented  in  the  following 
section . 


b.  Method  of  Thuras,  Jenkins,  and  O'Neil 

Thuras  et  al.  suggested  an  approximate  method  to  take  account 
of  ordinary  absorption  in  the  initial  growth  of  p^.  The  basis  for  their 
method  is  the  second  order  perturbation  solution  of  Eq.  2-1, 


(3kpi02x 

p2  ~ 2 
2p  c 
Ko  o 


The  spatial  growth  rate  (dp^dx)g  is,  from  Eq.  2-20, 


(»).• 


g 2p  c 
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(2-20) 


(2-21) 


Thuras  et  al.  reasoned  that  in  the  presence  of  absorption,  p^Q  should  be 
replaced  by  p^Qexp(-a^x)  where  is  the  absorption  coefficient  at  the 


fundamental  frequency.  From  this  reasoning  Eq.  2-21  becomes 
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(2-22) 


Thuras  et  al.  reasoned  that  opposing  the  spatial  growth  rate  represented 
by  Eq.  2-22  was  a spatial  decay  (dp^dxj^-a^p^  due  to  ordinary  absorption. 
Here  a 2 is  the  absorption  coefficient  at  the  frequency  of  the  second 
harmonic.  The  total  growth  rate  dp^dx  was  taken  to  be 
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The  solution  of  Eq.  2-23  which  satisfies  p2=0  at  x=0  is 
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(2-24) 


In  Chapter  3,  Section  B,  we  shall  make  use  of  Eq.  2-24  for  the 
expression  for  p2  required  by  Merklinger's  model.  It  is  interesting  to 
note  that  if  we  set  a2=4a^  (as  would  be  the  case  for  a viscous  fluid) 

Eq.  2-15  results. 


c.  Rudnick  Model 
24 

Rudnick  proposed  a model  for  describing  the  decay  of  sawtooth 
waves  in  a tube.  He  assumed  that  the  total  decay  rate  of  the  peak 
pressure  is  the  sum  of  the  decay  rate  due  to  small-signal  absorption 
and  the  decay  rate  due  to  nonlinear  effects  (found  from  Eq.  2-5)-  Thus 
his  model  is  similar  to  Thuras,  Jenkins,  and  O'Neil's  for  p2«  Rudnick 
obtained 


a P - 

o p 


P2 
2 P 


*poCo 


(2-25) 


where  a is  a "zero  order"  coefficient  for  the  tube  wall  attenuation  of 
o 
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Hie  solution  of  Eq.  2-25  consistent  with  F^=P^(0)  at  x=0  is 


the  wave 
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A modified  version  of  Eq.  2-25  is  discussed  in  Chapter  5 as  the  theoretical 
model  for  the  decay  of  the  fundamental  component  of  a sawtooth  wave,  and 
also  in  Appendix  C for  the  spherical  wave  problem.  Equation  2-25  is  quite 
appealing  physically  and  also  yields  simple  solutions.  As  is  always  the 
case,  however,  the  telling  test  is  the  experiment. 

4.  Summary 

We  have  discussed  several  theoretical  approaches  that  might  be 
used  in  the  description  of  saturation  of  plane  waves  in  a tube.  The 
Burgers'  equation  method,  though  certainly  applicable,  suffers  from  a lack 
of  appropriate  analytic  solutions.  The  computational  method  of  Pestorius, 
though  appropriate,  is  not  efficient  for  computing  saturation  curves.  The 
ad  hoc  models  do  not  include  the  effect  of  dispersion  but  they  have  yielded 
analytic  solutions.  The  models  pertain  to  the  fundamental  component  p . 
Merklinger's  model  is  based  on  an  intensity  equation,  and  Rudnick's  model, 
on  an  amplitude  equation.  We  shall  use  the  ad  hoc  models  for  our  theo- 
retical models  for  p^.  The  telling  test  of  the  models  is,  of  course,  the 
experiment . 


CHAPTER  5 


I 


THEORETICAL  DEVELOPMENT 


A.  Introduction 

In  this  chapter  the  theoretical  models  used  to  describe  the 
approach  to  saturation  are  presented.  Attention  is  focused  on  the  funda- 
mental frequency  component  p^  of  the  wave,  though  saturation  is  the 
ultimate  fate  of  all  the  Fourier  components.  First,  weak-shock  theory  is 
used  to  obtain  formulas  for  saturation  and  the  approach  to  saturation. 

We  then  endeavor  to  include  the  effects  of  tube  wall  dissipation. 

Although  dispersion  should  be  expected  to  play  a role  in  determining  the 
harmonic  amplitudes,  we  have  found  no  means  to  include  it  analytically. 

We  do,  however,  make  use  of  the  algorithm  of  Pestorius  to  give  an  estimate 
of  the  effect  of  dispersion.  This  algorithm  is  not  used  to  compute  actual 
saturation  levels  because  of  the  computation  time  required. 


B.  Estimates  of  the  Approach  to  Saturation 
1.  Weak-Shock  Theory 

An  estimate  of  the  saturation  amplitude  may  be  found  from 
weak-shock  theory.  The  expression  for  the  fundamental  pressure  p^  in  the 
sawtooth  region  is,  from  Eq.  2-5, 
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The  saturation  amplitude  Pl8at*  found  by  letting  the  source  amplitude  p1Q 
become  large  (i.e.,  o»l),  is 
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Thiis  P^sa^.  is  predicted  to  decrease  with  both  frequency  and  distance.  The 
inverse  first  power  dependence  of  P^ga^.  on  both  frequency  and  distance  is 
found  only  in  cases  where  absorption  over  the  entire  waveform  is  unimpor- 
tant. It  is  generally  true,  however,  that  the  saturation  level  does 
decrease  with  increasing  frequency  and  distance. 


2.  Merklinger's  Model 

As  was  mentioned  in  Chapter  2,  Merklinger’s  model  may  be 
reinterpreted  so  that  it  applies  to  more  than  just  viscous  fluids.  To  do 
so,  we  make  use  of  Thuras,  Jenkins,  and  O'Neil's  expression  for  the  second 
harmonic  pressure  p2  (Eq.  2-24).  Substituting  Eq.  2-24  for  Eq.  2-15  and 
following  the  same  analysis  as  in  Chapter  2,  Section  3. a.,  we  find  the 
differential  equation  for  p^  to  be 
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The  solution  of  Eq.  3-3  subject  to  the  boundary  condition  p^(0)=p^Q 
simply 


If  the  small-signal  attenuation  is  mainly  due  to  tube  wall  effects, 
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0*2= y2a^.  Note  that  if  as  would  be  the  case  for  a viscous  fluid, 

Eq.  3-4  reduces  to  Merklinger's  result,  Eq.  2-19. 

What  about  the  region  of  validity  of  Eq.  3-3?  We  should  expect 
that  use  of  the  approximation  p^p1+p2  (Eq.  2-10)  and  the  perturbation 
result  for  p2  (Eq.  2-21)  would  limit  the  use  of  Eq.  3-^  to  weak  waves. 
Equation  3-4  seems  to  have  wider  validity,  however.  Let  us  examine,  for 
example,  the  saturation  level  found  from  Eq.  3-4.  The  result  is 
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Now  let  ai,a2->0.  The  weak-shock  theory  expression,  Eq.  3-2,  is  recovered. 
We  conclude  that  although  several  of  the  assumptions  leading  to  Eq.  3-4 
are  not  valid  for  strong  waves,  the  saturation  level  is  an  appropriate 
strong-wave  limit. 


3*  The  Rudnick  Model 

We  describe  here  a generalization  of  the  method  Rudnick  used  to 
obtain  a formula  for  the  peak  amplitude  of  a sawtooth  wave.  The  "Rudnick 
assumption"  is  that  the  decay  rate  of  the  fundamental  pressure  component 
dp^/dx  is  the  sum  of  the  decay  rate  due  to  absorption  (dp^/dx)^^  =-<*., 
and  the  decay  rate  (dp ,/dx)  due  to  finite -amplitude  effects.  The 


latter  decay  rate  is,  from  Eq.  3-1, 
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The  total  decay  rate  dp^/dx  Is  then 
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The  solution  of  Eq.  3-7  is 
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where  p^O)  is  the  fundamental  pressure  amplitude  at  x=0.  If  we  require 
that  Eq.  3-8  reduce  to  the  weak-shock  solution  (Eq.  3-1)  as  a^+0,  then 
p^(0)  must  have  the  value 

P]_(0)  = 2p10  . (3-9) 


This  may  seem  a rather  surprising  result,  since  we  have  previously  speci- 
fied that  the  source  amplitude  is  P^q>  not  2p^Q.  The  discrepancy  is  only 
apparent,  however,  and  not  real.  If  we  extrapolate  the  sawtooth  solution 
(Eq.  2-5)  back  to  x=0  we  obtain, 

p(0,t)  = 2p10  ^ i sin  nout  , (5-10) 

which  represents  a sawtooth  wave  of  amplitude  np^Q.  This  suggests  that 
for  distances  a>3  the  acoustic  signal  produced  by  a sine  wave  of  amplitude 
p10  at  x=0  is  the  same  as  that  produced  by  a sawtooth  wave  whose  funda- 
mental has  amplitude  2p^Q  at  x=0.  Clearly  thiB  will  be  the  case  only  if 
absorption  is  not  important  over  the  initial  propagation  distance  of  three 


25 


shock  formation  lengths.  We  shall  presently  investigate  the  conditions 
under  which  this  assumption  is  reasonable. 

The  equation  for  p^  may  be  found  by  combining  Eq.  3-8  and 
Eq.  3-9-  The  result  is 
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Again,  as  the  source  amplitude  p^  increases  (i.e.,  x-*0),  the  saturation 
amplitude  p is  obtained, 
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It  is  seen  that  the  saturation  amplitude  still  decreases  with  increasing 

frequency  and  distance,  though  not  in  the  same  inverse  relation  given  by 

the  weak-shock  formula,  Eq.  3-2. 

We  now  investigate  the  conditions  under  which  the  Rudnick  model 

should  be  valid.  In  general,  Bince  the  model  equation,  Eq.  3-7  > is 

partially  derived  from  the  sawtooth  solution,  an  inherent  assumption  is 

that  the  wave  be  strong  enough  to  form  a sawtooth.  But  when,  in  the 

presence  of  ordinary  absorption,  is  sawtooth  formation  likely?  (See 

Appendix  A for  a similar  discussion  for  the  spherical  wave  problem. ) The 

shock  and  sawtooth  formation  distances  in  the  absence  of  absorption  may  be 

easily  computed  for  a given  set  of  source  and  medium  parameters.  What  is 

needed  is  an  estimate  of  the  distance  x beyond  which  small-signal 

max 

dissipation  effects  are  more  important  than  nonlinear  dissipation  effects. 
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Blackstock  computed  the  distance  x at  which  the  decay  rates  due  to 

max  J 

nonlinear  effects  (found  from  Eq.  5-l)  and  small-signal  effects  are  equal. 

For  plane  waves  the  value  of  x is 

max 
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The  question  of  shock  and  sawtooth  formation  may  thus  be  answered  by 

_ /\  ^ 

analyzing  the  various  possible  relations  between  x,  x,  and  x .'"If  the 

’ max 

relation  between  x and  x is 
max 
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shock  formation  is  precluded.  The  high  frequency  components  damp  out 
before  a shock  can  form.  The  inequality 


x < x 
max 


(3-15) 


indicates  that  shock  formation  i6  possible.  Ordinary  dissipation  domi- 
nates, however,  before  a sawtooth  can  form.  Sawtooth  formation  is  likely 
if  the  relation  is 


max 


(3-16) 


We  see  that  the  Rudnick  model  should  be  valid  if  Eq.  3-l6  is  satisfied. 
That  is,  Eq.  5-11  should  be  valid  if  x is  less  than  l/a. 


CHAPTER  4 


EXPERIMENTAL  APPARATUS  AND  PROCEDURES 


A.  Introduction 

In  this  chapter  the  equipment  and  procedures  used  in  the 
experiments  described  in  Chapter  5 are  discussed.  Individual  equipment 
listings  are  given  in  Section  B along  with  a general  schematic  of  the 
experimental  apparatus.  In  Section  C the  experimental  procedures  axe 
discussed. 

B.  Experimental  Apparatus 

Figure  4-1  shows  a generalized  schematic  of  the  experimental 
apparatus.  The  system  has  been  used  in  several  previous  experiments  (see, 
for  example,  Refs.  8 and  15).  The  microphone  nearest  the  driver,  called 
a source  or  monitor  microphone,  is  used  to  measure  the  source  waveform. 
This  microphone  also  provides  an  input  to  an  automatic  gain  control  (AGC) 
circuit,  which  is  used  to  stabilize  the  source  SPL.  The  second  micro- 
phone is  the  receiving  microphone,  which  is  variable  in  5-7  m steps 
relative  to  the  source  microphone. 

The  experimental  apparatus  is  divided  into  three  systems: 
a transmitting  system,  the  plane  wave  tube,  and  the  receiving  system. 

The  specific  pieces  of  equipment  used  in  each  of  these  systems  are  listed 
in  the  following  three  sections. 
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FIGURE  4-1 

BLOCK  DIAGRAM  OF  THE  EXPERIMENTAL  APPARATUS 
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1.  The  Transmitting  System 

The  following  pieces  of  equipment  were  used  in  the  transmitting 

system. 

(1)  Oscillator.  B&K  type  1022  or  H-P  3300A  with  3302  A 
trigger/ phase  lock  plug-in. 

(2)  Counter.  H-P  5300  B. 

(3)  Gate.  GR  1496. 

(4)  Amplifier.  Altec  model  250  B power  amplifier. 

(5)  Oscilloscope.  Tektronix  532. 

(6)  Driver.  JBL  model  375-H.  This  driver  is  commercially 
available  with  either  a phenolic  or  aluminum  diaphragm. 
Because  of  its  much  better  high  frequency  response,  the 
driver  with  the  aluminum  diaphragm  was  chosen  for  this 
experiment.  A typical  frequency  response  of  one  of  the 
drivers  is  shown  in  Fig.  4-2.  The  response  wa6  taken 


FIGURE  4-2 

TYPICAL  FREQUENCY  RESPONSE 
OF  JBL  375-H  DRIVER 

(CONSTANT  VOLTAGE  INPUT) 
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with  a constant  voltage  input  to  the  driver.  For 
this  measurement  the  driver  was  coupled  to  the  plane 
wave  tube  shown  in  Fig.  4-1.  What  is  plotted  in 
Fig.  4-2  is  the  relative  SPL  measured  at  the  monitor  micro- 
phone, which  is  approximately  13  cm  from  the  driver  throat. 
The  response  of  the  driver  is  flat  within  a 3 dB  range  from 
approximately  300  Hz  to  85OO  Hz.  The  maximum  input  power 
to  the  driver  is  specified  as  30  Wrms.  It  was  found, 
however,  that  driver  failure  occurred  much  less  frequently 
when  the  input  power  was  limited  to  25  Wrms.  At  25  Wrms 
input  the  maximum  SPL  at  the  monitor  microphone  was  157  dB 
at  500  Hz.  To  attain  higher  levels,  provision  was  made  for 
parallel  operation  of  four  375_H  drivers.  The  maximum  SPL 
at  the  first  microphone  was  then  163  dB  at  500  Hz. 

(7)  Laboratory  constructed  predistortion  network  (not  shown  in 
Fig.  4-1).  This  network  was  used  to  predistort  the  elec- 
trical input  to  the  driver  to  compensate  for  driver 
distortion  and  for  finite-amplitude  propagation  distortion 
over  the  13  cm  path  to  the  monitor  microphone.  The  network 
is  discussed  in  detail  in  Chapter  5« 

2.  Plane  Wave  Tube 

The  tube  used  in  this  experiment  was  originally  constructed  by 

McKittrick;  see  Ref.  l6.  It  consists  of  eight  3.7  m sections  of  aluminum 

pipe  Joined  by  microphone-holding  flanges.  The  pipe  has  a 2 m long  fiber- 
o 

glass  termination.  A typical  flange  assembly  is  shown  in  Fig.  4-3,  taken 
from  Ref.  15.  The  l/4  in.  microphone  fits  through  the  white  teflon  collar 


FIGURE  4-3 

EXPLODED  VIEW  OF  FLANGE  ASSEMBLY  WITH 
MEASURING  MICROPHONE  AND  PLUG 

(THE  PROTECTIVE  GRID  SHOWN  ON  THE  MICROPHONE 
WAS  REMOVED  BEFORE  THE  MICROPHONE  WAS 
INSERTED  IN  THE  HOLDER  AND  FLANGE) 


Taken  from  Ref.  15 
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to  mount  flush  with  the  inside  surface  of  the  tube  wall.  When  the 
microphone  was  not  used,  an  aluminum  plug  (shown  directly  below  the  micro- 
phone in  Fig.  4-3)  was  installed.  The  termination  was  designed  by 

8 23 

Williams,  on  the  basis  of  a paper  by  Burns.  The  first  part  of  the 

termination  was  a tapered  section  about  1 m long.  The  second  part,  also 
1 m long,  completely  filled  the  pipe  cross  section.  The  reflection 
coefficients  of  the  termination  were  measured  by  a pulse  method;  the 
results  are  given  in  Chapter  5* 

3.  The  Receiving  System 

The  following  pieces  of  equipment  were  used  in  the  receiving 

system. 

(1)  Microphone.  B&K  type  4136,  l/4  in.  The  pressure  response 
of  this  microphone  is  flat  to  within  ±0.5  dB  from  50  Hz  to 
70  kHz.  When  flush  mounted  in  the  tube  wall,  however,  the 
microphone  acts  as  a baffled  circular  piston  receiver. 
Because  the  sound  waves  are  incident  at  90°  to  the  piston, 
the  3 dB  down  frequency  response  of  the  microphone  i6 
limited  to  30  Hz  to  28  kHz.  All  other  components  in  the 
receiving  system  have  a 3 dB  down  response  of  at  least 

20  Hz  to  200  kHz.  The  microphone  is  therefore  the  limiting 
element  in  the  receiving  system. 

(2)  Preamplifier.  B&K  type  2619. 

(3)  Preamplifier  and  Power  Supply.  B&K  type  2803  dual  channel. 

(4)  Microphone  Cables.  B&K  type  AD0029  30  m cables. 

(5)  Microphone  Calibrator  (not  shown  in  Fig.  4-l).  B&K  type 
4220  pistonphone.  The  output  of  the  calibrator  is  a 240  Hz 
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sinusoid  with  a SPL  of  123.8  *0.2  dB. 

(6)  Spectrum  Analyzer.  H-P  model  3580A. 

(7)  Oscilloscope  Camera  (not  shown  in  Pig.  4-1). 

H-P  model  197- A. 

(8)  Voltmeter.  H-P  model  400EL. 

(9)  Oscilloscope.  Tektronix  model  545B. 

(10)  Oscilloscope  Camera  (not  shown  in  Fig.  4-1).  Tektronix 
model  125 • 

(11)  Preamplifier.  This  preamplifier  was  constructed  for  use 
as  a wideband  amplifier  with  negligible  phase  distortion. 
The  3 dB  down  frequency  response  of  this  amplifier  is  7 Hz 
to  700  kHz.  There  is  negligible  phase  distortion  above 
100  Hz.  The  gain  of  the  amplifier  iB  20  ±0.1  dB  from 

500  Hz  to  200  kHz. 

C.  Experimental  Procedures 

Three  experiments  were  conducted  in  this  investigation.  The 
procedures  are  discussed  below. 

1.  Measurement  of  Small-Signal  Attenuation  Coefficients 

The  procedure  here  was  to  generate  a low  level  (SPL  = 115  to 
120  dB)  sinusoid  at  the  source  and  observe  the  decay  of  the  signal  with 
distance.  The  attenuation  coefficient  at  a particular  frequency  was  then 
obtained  from  a least  squares  fit  of  the  plot  of  measured  loss  versus 


distance 
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2.  Measurement  of  Pipe  Termination  Reflectivity 

The  pulse  method  was  used  to  measure  the  reflectivity  of  the 
pipe  termination.  The  receiving  microphone  was  placed  approximately  1.7  m 
from  the  beginning  of  the  pipe  termination.  The  direct  (pfl)  and  reflected 
(Pr)  pulse  amplitudes  were  observed  on  the  oscilloscope  and  their  respec- 
tive values  recorded;  the  amplitude  of  the  reflected  pulse  was  corrected 
for  the  round  trip  loss  caused  by  tube  wall  attenuation.  The  reflection 
coefficient  R was  determined  from  the  expression 

r 2c  *i 


R = 20  log10 


(4-1) 


3»  Measurements  of  Received  Acoustic  Signal  versus  Propagation 

Distance  and  Source  Level 

With  the  available  apparatus  we  could,  in  principle,  vary  p1Q, 
u>,  and  x.  The  most  convenient  way  of  testing  the  various  theoretical 
predictions  of  saturation  was  to  fix  x and  <n  and  vary  p^Q.  In  this  way  a 
family  of  amplitude  response  curves  was  obtained.  For  the  purpose  of 
investigating  certain  other  effects,  such  as  dispersion,  it  was  preferable 
to  hold  p^Q  and  co  fixed  and  to  vary  x.  In  any  case  the  procedure  was  to 
measure  the  received  signal  as  a function  of  x or  P^q*  Each  measurement 
consisted  of  recording  an  oscillogram  and  spectrum  of  the  received  signal. 

4.  Measurement  Accuracy 

The  accuracy  associated  with  the  measurement  of  SPL,  f,  and  x 
was  estimated  from  the  various  uncertainties  in  the  response  of  the  trans- 
mitting and  receiving  system.  The  SPL  measurements  are  believed  accurate 


to  ±0.5  dB.  The  frequency  stability  of  the  B&K  oscillator  was  observed  to 
be  ±0.2%.  Distance  measurements  are  accurate  to  ±1$. 
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CHAPTER  5 

EXPERIMENTAL  AND  THEORETICAL  RESULTS 


A.  Introduction 

In  this  chapter  the  experimental  results  obtained  are  given  and 
compared  with  the  theoretical  predictions  found  in  Chapter  J.  These 
results  are  divided  into  two  sections,  the  first  dealing  with  the  salient 
acoustical  properties  of  the  pipe,  the  second  with  the  high  intensity 
experiment.  All  computations  were  made  for  the  measured  temperature, 
relative  humidity,  and  ambient  pressure  of  20°C,  50$,  and  760  mm  Hg., 
respectively. 

Some  of  the  data  given  here  have  been  previously  reported. 
Figures  5-5,  5-8,  and  5-10  were  reported  in  Ref.  26.  Figure  5-5  was 
reported  in  Ref.  27. 

B.  Acoustical  Properties  of  the  Pipe 

1.  Measurement  of  Attenuation  Coefficients 

Figure  5-1  shows  the  measured  frequency  dependence  of  the 
attenuation  coefficient  for  the  pipe.  The  points  are  the  experimental 
data  and  the  dashed  line  is  a least  squares  fit  to  the  data.  The  solid 
line  is  the  attenuation  given  by  the  Kirchhoff  formula  (Eq.  2-6).  The 
offset  of  the  two  lines  is  about  10$.  This  is  similar  to  the  findings  of 
Pernet  and  Payne, ^ (8$),  and  Pestorius,®  (10$). 
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2.  Measurements  of  the  Termination  Reflectivity 

As  was  mentioned  in  Chapter  4,  the  reflection  coefficients  of 
the  termination  were  measured  by  the  pulse  method.  A plot  of  R,  deter- 
mined from  Eq.  4-1,  versus  frequency  is  shown  in  Fig.  5-2.  The  reflection 
coefficients  are  typical  for  such  a termination.  We  conclude  that  for  the 
frequency  range  shown  the  tube  is  properly  terminated. 

C.  Measurements  of  Initially  Sinusoidal  Waves 
1.  Predistortion 

In  Chapter  4 it  was  noted  that  at  high  amplitudes  the  acoustic 
signal  at  the  monitor  microphone  would,  in  some  cases,  be  distorted.  This 
distortion  is  a combination  of  distortion  in  the  horn  driver  itself  and 
propagation  distortion  in  the  13  cm  separating  the  horn  driver  and  the 
monitor  microphone.  In  all  cases  the  distortion  was  primarily  second 
harmonic.  One  would  expect  that  the  electrical  input  to  the  driver  could 
be  predistorted  in  such  a way  as  to  minimize  the  distortion  of  the  acoustic 
wave  arriving  at  the  monitor  microphone.  An  electrical  network  designed 
to  produce  the  predistorted  signal  is  shown  in  Fig.  5-3.  The  oscillator 
produces  a sine  wave  at  the  second  harmonic  of  the  desired  operating  fre- 
quency. Hiis  signal  is  fed  into  one  input  of  a summing  amplifier  and  also 
into  a divide-by-two  counter.  The  counter  drives  the  phase-locked  oscil- 
lator. Hie  phase-locked  oscillator  produces  a sine  wave  at  the  operating 
frequency  and  serves  as  the  second  input  to  the  summing  amplifier.  The 
output  of  the  summing  amplifier  is  then  the  sum  of  signals  at  the  funda- 
mental and  second  harmonic  frequencies.  The  phase  of  the  second  harmonic 
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BLOCK  DIAGRAM  OF  PREDISTORTION  NETWORK 
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component  may  be  varied  ±180°  with  respect  to  the  fundamental  by  means  of 
a phase  control  on  the  phase-locked  oscillator. 

Typical  results  using  the  predistortion  network  are  shown  in 
Fig.  5-4.  It  is  evident  that  although  the  electrical  input  to  the  driver 
is  a sinusoid  the  acoustic  signal  is  distorted.  Spectral  analysis  of  the 
acoustic  signal  shows  the  presence  of  a strong  second  harmonic  component 
15  dB  below  the  fundamental.  The  second  column  of  Fig.  5-4  shows  the 
results  of  using  predistortion.  When  the  phase  and  amplitude  of  the 
second  harmonic  fed  to  the  driver  are  varied  intil  the  second  harmonic  is 
minimized  at  the  monitor  microphone,  a reduction  of  22.5  dB  in  the  acous- 
tic second  harmonic  is  achieved. 

The  main  limitation  of  using  predistortion  is  the  bandwidth  of 
the  electroacoustic  transducer.  In  our  case  bandwidth  is  no  limitation 
because  the  upper  5 dB  down  frequency  of  the  driver  is  8.5  kHz  and  our 
highest  fundamental  frequency  is  4 kHz.  With  suitable  counters  the  method 
could  be  used  to  suppress  any  harmonic  within  the  transducer  bandwidth. 

2.  Amplitude  Response  Curves 

The  phenomenon  of  acoustic  saturation  can  most  easily  be  seen 
from  an  amplitude  response  curve  taken  at  a fixed  distance.  Two  such 
curves  are  shown  in  Fig.  5-5.  The  received  SPL  at  the  fundamental  fre- 
quency is  plotted  versus  source  level  for  two  propagation  distances  at 
several  source  frequencies.  The  dashed  lines  are  predictions  based  on 
linear  theory,  and  the  solid  curves  (valid  for  c>5)  are  the  predictions 
based  on  the  Rudnick  model  (Eq.  3-12).  Saturation  is  well  evidenced  by 
the  bending  over  of  the  response  curves  at  high  source  levels.  For 
example,  when  f = 3.57  kHz  and  x = 25.8  m,  the  extra  attenuation  is  about 
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1^.1  dB.  At  this  point  96$  of  the  power  at  the  fundamental  frequency  has 
been  lost  because  of  nonlinear  effects! 

At  the  14.8  m distance  the  theoretical  curves  are  in  good 
agreement  with  the  experimental  data.  At  25.8  m,  however,  the  theoretical 
curves  fall  slightly  below  the  data,  the  agreement  becoming  worse  as  the 
frequency  increases.  The  reason  for  this  small  discrepancy  is  discussed 
in  Chapter  6. 

We  have  thus  far  discussed  saturation  in  terms  of  the 
fundamental  frequency  only.  Saturation  is,  however,  the  ultimate  fate  of 
each  Fourier  component.  Figure  5-6  shows  an  amplitude  response  curve  for 
the  first  three  harmonics  of  an  Initially  sinusoidal  wave.  Hie  solid 
curves  are  best  fitB  to  the  data.  A theoretical  curve  was  not  Included 
for  the  fundamental  component  because  these  data  were  compared  with  theory 
in  Fig.  5-5.  We  have  no  theory  for  the  higher  harmonics  comparable  to 
that  for  the  fundamental.  We  may,  however,  see  from  the  weak-shock 
formula  (Eq.  2-5)  that  the  higher  harmonic  components  do  Indeed  saturate. 
This  result  is  borne  out  by  Fig.  5-6.  The  approach  to  saturation  of  each 
harmonic  is  evident. 

Hie  other  models  discussed  in  Chapter  3 were  also  compared  to 
data.  A typical  result  is  shown  in  Fig.  5-7.  The  equation  based  on  weak- 
shock  theory  (Eq.  3-1)  yields  a curve  about  5 dB  above  the  data.  Ordinary 
absorption  is  Indeed  important  here.  The  small-signal  loss  for  this  fre- 
quency and  distance  is  about  7.5  dB.  The  equation  based  on  Merkllnger's 
model  provides  a good  fit  to  the  data  in  the  linear  region.  In  the 
plateau  region,  however,  it  yields  a curve  about  1.3  dB  above  the  data. 

Hie  equation  based  on  the  Rudnick  model  fits  the  data  to  within  0.1  dB. 


A. 


FIGURE  5-6 

AMPLITUDE  RESPONSE  CURVES  FOR 
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AN  INITIALLY  SINUSOIDAL  WAVE 
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COMPARISON  OF  THE  VARIOUS  THEORETICAL 
MODELS  WITH  DATA 
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In  all  cases  the  Rudnlck  model  equation  provided  a best  estimate  of  the 
approach  to  saturation.  For  this  reason  the  Rudnlck  model  was  used  as  our 
theoretical  model  for  the  rest  of  this  work.  The  fact  that  the  Rudnlck 
model  is  limited  to  waves  strong  enough  to  form  a sawtooth  Is  not  critical 
for  our  work  because  we  are  interested  in  saturation  and  the  approach  to 
saturation . 

3*  Time  Waveforms  of  Received  Signals 

The  received  waveforms  as  a function  of  source  SPL  are  also  of 
interest.  (See  Fig.  5-8* ) The  waveforms  at  high  levels  illustrate  the 
strong  distortion  that  is  responsible  for  the  extra  attenuation.  Satura- 
tion is  evident  by  inspection  of  the  last  two  waveforms.  The  amplitudes 
are  almost  the  same  despite  a 3.2  dB  difference  in  the  source  SPL.  It 
seems  likely  that  the  minor  "raggedness"  in  the  lower  three  curves  is  due 
to  scattering  from  small  irregularities  in  the  pipe  near  the  microphone. 
These  irregularities,  probably  associated  with  the  flange  Junction  and/or 
the  microphone  mounting  hole,  would  be  expected  to  scatter  only  very  high 
frequencies.  Such  an  expectation  is  consistent  vlth  the  fact  that  the 
raggedness  appears  only  in  waveforms  that  contain  shocks.  Note  the  asym- 
metry of  the  waveforms.  It  begins  to  show  up  at  142  dB  and  is  very 
apparent  at  148  dB.  The  asymmetry  is  caused  by  dispersion.  The  asymmetry 
is  lessened,  however,  as  the  source  amplitude  is  increased.  See  the  last 
two  waveforms.  Nonlinear  effects  have  become  so  strong  that  "steepening” 
of  the  shock  tends  to  overcome  the  rounding  of  the  shock  peak. 

Modified  weak-shock  theory  is  used  in  Chapter  6 to  investigate 
the  effect  of  dispersion.  It  is  therefore  appropriate  here  to  compare  our 
experimental  waveforms  vlth  those  computed  using  Festorlus'  algorithm. 


SOURCE  SPL 


130.0  dB 


136  0 dB 


OSCILLOSCOPE 

GAIN 

6 dB 


6 dB 


0 dB 


0 dB 


0 dB 


0 dB 


FIGURE  5*8 

OSCILLOGRAMS  OF  RECEIVED  SIGNAL  FOR 
VARIOUS  SOURCE  LEVELS  ( f-2  kHz,  »-14.8  m) 


49 


Figure  5-9  shows  such  a comparison.  The  top  waveform  in  each  column  is 
the  same  because  the  measured  source  waveform  was  used  as  the  source 
waveform  for  Pestorius'  algorithm.  The  absorption  coefficient  used  in 
the  computer  algorithm  is  the  measured  coefficient  from  Fig.  5-1.  We  note 
that  the  agreement  between  experiment  and  theory  is  quite  good  in  all 
respects.  The  progressive  distortion  of  the  wave  with  propagation  dis- 
tance, which  is  responsible  for  saturation,  is  also  evident  in  Fig.  5-9. 

At  3*7  m the  waveform  is  a sawtooth,  though  it  is  slightly  asynsnetrlc. 

As  the  distance  increases  the  asymmetry  grows;  the  shock  peak  becomes  more 
rounded  while  the  trough  remains  sharp.  The  reason  is  that  as  the  dis- 
tance increases  the  wave  amplitude  is  steadily  diminished.  Nonlinear 
effects  become  too  weak  to  overcome  the  rounding  of  the  shock  peak  caused 
by  dispersion. 

4.  Propagation  Curves 

Figure  5-10  shows  propagation  curves  for  the  fundamental 

component  for  several  different  source  frequencies.  Hie  solid  curves  are 

predictions  based  on  the  Rudnlck  model  (Eq.  3-ll)>  whereas  the  dashed 

curves  are  those  based  on  weak-shock  theory  (Eq.  3-l) • Each  of  these 

curves  stops  at  the  left  at  the  distance  x=x  because  the  predictions  see 

not  valid  at  lesser  distances . The  dotted  line  at  the  top  is  a small- 

signal  prediction  for  the  500  Hz  case.  Hie  point  labeled  x is  the 

max 

distance  beyond  which  ordinary  absorption  is  more  important  than  shock 
dissipation.  Beyond  x weak-shock  theory  is  clearly  not  valid.  We 
note  here  that  for  the  data  shown  in  Fig.  5*10,  the  weak-shock  formula 
overestimates  the  data  by  about  5 dB  at  x°xm^v . Hie  Rudnlck  equation 
appears  to  be  a great  improvement  over  the  weak- shock  formula. 
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CHAPTER  6 


CONCLUSIONS 


A.  Introduction 

In  this  chapter  the  results  are  summarized  and  conclusions  are 
drawn.  First  the  role  of  dispersion  is  discussed.  Then  the  success  of 
each  of  the  various  models  used  to  predict  saturation  and  the  approach  to 
saturation  is  examined.  Finally,  general  conclusions  about  the  research 
are  given. 

B.  The  Role  of  Dispersion 

In  the  Rudnick  model  introduced  in  Chapter  2,  dispersion  is  not 

taken  into  account.  As  shown  in  Chapter  2,  dispersion  is  responsible  for 

the  asymmetric  waveforms  observed  in  the  pipe.  The  shock  peak  is  rounded 

while  the  trough  remains  sharp.  Since  shock  dissipation  depends  on  the 

actual  magnitude  of  the  pressure  discontinuity  at  the  shock,  the  rounding 

of  the  shock  peak  reduces  shock  dissipation.  In  other  words,  dispersion 

o 

inhibits  the  attenuation  due  to  shocks.  By  ignoring  dispersion  in  our 
theoretical  model,  Eq.  3-11,  we  have  overemphasized  shock  dissipation  and 
have  therefore  arrived  at  predicted  levels  that  are  too  low.  Furthermore, 
since  the  phase  shift  caused  by  dispersion  depends  on  the  parameter  ctx, 
we  should  expect  the  discrepancy  between  prediction  and  experiment  to 
increase  with  both  frequency  and  propagation  distance.  This  expectation 
is  qualitatively  consistent  with  the  data  and  curves  in  Fig.  5-5  and 
Fig.  5-10. 
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We  have  argued  here  that  dispersion  could,  qualitatively,  be 
responsible  for  the  slight  discrepancies  observed  between  data  and  the 
Rudnick  model  at  the  greater  distances  and  higher  frequencies.  Pestorius' 
algorithm  was  used  to  find  the  quantitative  effect  of  dispersion  on  the 
received  signal.  A check  on  the  computer  algorithm  was  made  by  comparing 
the  predicted  propagation  curve  with  experimental  data  from  the  tube. 

(See  Fig.  6-1.)  The  computation  was  made  for  the  case  of  a 3*5  kHz  plane 
wave  with  a source  SPL  of  157  dB;  the  experimental  value  of  the  absorption 
coefficient  used  was  a = 0.078  Np/m.  The  agreement  with  the  data  is 
excellent;  the  maximum  deviation  is  0.5  dB.  To  estimate  the  effect  of 
dispersion  on  the  signal,  an  additional  computer  run  was  made  ignoring 
dispersion  effects.  The  propagation  curve  obtained  from  this  run  is  not 
shown  in  Fig.  6-1  because  it  differs  at  most  by  0.1  dB  from  the  curve 
shown  there.  We  conclude  that  for  this  case  dispersion  has  a negligible 
effect  on  the  amplitude  of  the  fundamental  component.  Since  the  given 
conditions  represent  the  highest  frequency  and  greatest  propagation  dis- 
tance in  the  experiment,  dispersion  alone  cannot  account  for  the  small 
discrepancies  between  the  Rudnick  model  and  the  measured  data. 

The  higher  harmonic  components  computed  in  the  two  runs  are 
compared  in  Table  6-1.  The  harmonic  amplitudes  with  dispersion  included 
(Bq)  are  compared  to  those  without  dispersion  (B'n).  The  predicted  funda- 
mental component  of  the  wave  is  indeed  higher  when  dispersion  effects  sure 
included;  however,  the  difference  is  negligible  for  the  given  conditions. 
Dispersion  does,  however,  seem  to  become  more  important  for  the  higher 
harmonics.  This  finding  is  consistent  with  the  observation  that  the 
shape  of  the  waveform  in  the  neighborhood  of  each  shock  is  strongly 
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FIGURE  6-1 

COMPARISON  OF  MEASURED  PROPAGATION 
DATA  WITH  PESTORIUS’  ALGORITHM 
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affected  by  dispersion.  The  shock  shape  is  largely  determined  by  the 
higher  harmonic  content. 

TABLE  6-1 

COMPARISON  OF  COMPUTED  HARMONIC  AMPLITUDES  WITH  (B  ) 

AND  WITHOUT  (B'  ) DISPERSION  CORRECTIONS  n 

Harmonic  Number  20  log^Q  [®n/®'n] 

1 0.1 
2 -0.01 

3 -0.2 

4-  -0.3 

5 -0.5 

6 -0.6 

C.  Comparison  of  The  Theoretical  Models 
1.  Modified  Weak -Shock  Theory 

Pestorius'  computer  algorithm  was  found  to  give  excellent 
agreement  with  the  data  in  all  the  cases  in  which  it  was  used.  It  has 
two  main  attractive  features.  First,  it  includes  the  elements  of  weak- 
shock  theory,  small-signal  absorption,  and  dispersion.  Second,  the  source 
waveform  that  may  be  considered  is  arbitrary.  The  fact  that  the  algorithm 
is  not  an  analytic  solution  made  it  poorly  adapted  to  our  use,  however. 

The  computation  of  an  amplitude  response  curve,  for  example,  would  have 
required  a separate  computer  run  for  each  increment  in  source  level. 


56 

2.  Modified  Merklinger  Model 

The  model  suggested  by  Merklinger  to  describe  the  fundamental 
component  of  an  initially  sinusoidal  plane  wave  in  a viscous  fluid  was 
generalized  so  that  it  could  be  applied  to  any  type  of  absorbing  fluid. 

The  basis  of  the  model  is  an  intensity  equation.  The  primary  advantage  of 
this  model  is  that  it  may  be  used  over  the  entire  range  of  source  level. 

In  all  cases  near  saturation,  however,  received  levels  predicted  using 
this  model  were  higher  than  those  actually  measured.  Typical  predicted 
results  are  shown  in  comparison  to  experimental  data  in  Fig.  5-7. 

5*  The  Rudnick  Model 

This  model  was  originally  used  by  Rudnick  to  describe  the  peak 
amplitude  of  a sawtooth  wave  in  the  presence  of  tube  wall  absorption. 

The  basis  of  the  model  is  an  amplitude  equation.  We  adapted  Rudnick' s 
approach  to  obtain  an  equation  for  the  decay  of  the  fundamental  component 
p^.  This  equation  was  solved,  and  predictions  were  found  to  be  in  good 
agreement  with  the  experimental  data.  Small  discrepancies  were  found, 
however.  These  discrepancies  increased  with  increasing  frequency  and 
distance.  Dispersion,  which  was  originally  thought  to  be  the  cause  of  the 
small  discrepancies,  was  found  to  be  insignificant  for  our  experimental 
conditions. 

D.  Summary  and  Conclusions 

The  primary  result  of  this  research  is  the  experimental  proof 
of  acoustic  saturation  of  plane  waves  in  air.  The  tendency  toward  satura- 
tion was  found  to  increase  with  source  frequency  and/or  distance.  In  one 
case  it  was  found  that  nonlinear  effects  caused  96#  of  the  power  at  the 
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fundamental  frequency  to  be  lost  over  a propagation  path  of  26  m.  The 
transfer  of  energy  to  higher  frequency  components,  which  is  responsible 
for  acoustic  saturation,  was  well  evidenced  by  the  formation  of  sawtooth 
waves.  Dispersion  radically  alters  the  time  waveform  of  the  received 
signal  (as  discussed  in  Section  B of  this  chapter) . However,  calculations 
based  on  modified  weak  shock  theory  show  that  for  our  experimental  con- 
ditions dispersion  has  negligible  effect  on  the  amplitude  of  the  fundamen- 
tal component. 

Another  important  result  of  this  research  is  the  finding  that 
existing  theoretical  models,  adapted  to  our  experimental  conditions,  gave 
a good  quantitative  explanation  of  saturation  and  the  approach  to  satura- 
tion. Though  weak-shock  theory  generally  indicated  the  development  of 
saturation,  it  was  necessary  to  take  explicit  account  of  ordinary  absorp- 
tion. We  adapted  the  methods  of  Rudnick  and  Merklinger  to  obtain 
expressions  for  the  fundamental  component  p^.  In  all  cases  the  Rudnick 
equation  was  in  better  agreement  with  data. 

Some  experimental  and  theoretical  results  for  spherical  waves 
are  presented  in  the  appendices. 
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In  some  cases  it  was  possible  to  generalize  the  theoretical 

work  done  for  plane  waves  to  the  case  of  spherically  spreading  waves . 

These  results  are  included  here  as  appendices.  In  Chapter  3,  for  example, 

the  question  of  whether  ordinary  dissipation  can  prevent  the  formation 

of  shocks  or  of  a fully  formed  sawtooth  was  answered  by  comparing  the 

critical  distances  x,  x,  and  x In  Appendix  A we  discuss  a method 

max 

for  assessing  the  importance  of  nonlinear  effects  on  spherical  waves 

based  on  a comparison  of  the  critical  ranges  r,  ?,  and  r . This 

max 

method  leads  to  a classification  of  nonlinear  effects  on  spherical 
waves  as  weak,  moderate,  or  strong.  A perturbation  solution  of  Burgers' 
equation  for  weak  spherical  waves  is  presented  in  Appendix  B.  An  exten- 
sion of  the  Rudnick  model  to  spherical  waves  i6  presented  in  Appendix  C. 
The  model  is  appropriate  for  strong  waves. 

The  weak  and  strong-wave  solutions  were  compared  with  data  from 
freefield  propagation  experiments.  In  both  the  weak  and  strong  wave 
experiments  the  measured  quantities  were  p^  and  p^  versus  range.  In  the 
weak-wave  experiment  the  source  was  an  acoustic  array;  for  the  strong- 
wave  experiments  an  acoustic  siren  was  used. 

We  note  here  that  the  material  presented  in  Appendix  A was 


included  in  Ref.  27. 
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APPENDIX  A 

ESTIMATES  OF  THE  IMPORTANCE  OF  NONLINEARITY 
ON  THE  PROPAGATION  OF  SPHERICAL  WAVES 


Here  we  outline  a method  for  estimating  the  importance  of 

nonlinear  effects  on  the  propagation  of  spherically  spreading  sound 

waves.  The  analysis  is  based  on  a comparison  of  the  three  critical 

range  parameters  r,  r,  and  r corresponding,  respectively,  to  shock 

max 

formation,  sawtooth  formation,  and  ’’old  age."  The  analysis  we  present 

here  is  somewhat  similar  to  that  of  Shooter,  Muir,  and  Blackstock.^  The 

difference  between  our  approach  and  that  of  Shooter  et  al . is  the 

particular  graphical  representation  of  the  resulting  equations. 

The  question  of  when  shock  or  sawtooth  formation  is  possible 

for  a plane  wave  was  dealt  with  in  Chapter  5-  The  analysis  for  spherical 

waves  is  precisely  the  same  if  we  replace  x,  x,  and  x with  r,  r,  and 

max 

rmax'  resPectlvely-  To  specific  formulas  for  r,  r,  and  r^,  we 

first  consider  the  case  of  an  omnidirectional  spherical  wave.  That  is, 
we  let  the  boundary  condition  be 

p = p1Q  sin  cot  , at  r = rQ  . ( A-l) 
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We  assume  here  that  kro»l.  The  critical  ranges  are^ 


1/pekr 


r = r e 
o 


3/p€kr 


r = r e 
o 


(A-2) 

(A-3) 


r 

max 


(A-4) 


To  simplify  the  assessment  of  the  importance  of  nonlinearity, 
we  now  consider  a graphical  display  of  Eqs.  A-2  through  A-1*.  A direct 
approach  would  be  to  plot  the  three  critical  distances  themselves.  If 
Eqs.  A-2  and  A-3  are  divided  by  rQ,  the  normalized  quantities  r/rQ  and 
*/T0  may  be  plotted  versus  the  single  parameter  0€kro<  The  equation  for 
rmax/ro>  however,  depends  on  two  parameters,  pekr^  and  Oro.  A family  of 
curves  would  therefore  have  to  be  plotted  to  represent  rmax/ rD • An 
alternative  approach  is  as  follows : a reference  point  roughly  marking 

the  threshold  of  importance  of  nonlinearity  is  *=rmax»  far  more  serious 
nonlinear  effects  are  indicated  by  r=rTnny.  These  two  equalities  define 
the  following  relations,  respectively; 

-l/p€kr 

? = rmax  "*  ^o  “ P€kro  6 ° ' (A“5) 


and 


max 


Uar 


pckro  e 


-3/p€krr 


(A-6) 


Only  two  dimensionless  groups  appear  in  Eqs.  A-5  and  A-6,  0ekr  and  Or  . 

o o 

This  suggests  representing  these  two  equations  by  plotting,  for  example. 
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pekr^  versus  ctr^.  It  is  easily  seen  that  for  a given  medium  the  quantity 
P€krQ  depends  only  on  the  source  level  SL  (farfield  SPL  extrapolated  to 
1 m)  and  f.  A convenient  measure  of  P€krQ  is  the  "source -frequency  level" 
SFL  (dB  re  0.0002  pbar-kHz-m)  defined  as  follows: 

SFL  = SL  + 20  log10  , (A-7) 

where  f^^  is  the  frequency  in  kilohertz.  The  SFL  as  we  have  defined  it 

here  is  what  Merklinger,  Mellen,  and  Moffet2^  call  the  "scaled  source 

level"  SL*.  The  relation  between  SFL  and  Bekr  for  air  at  20°C  is 

o 

SFL  = 167.2  + 20  log1Q  6€krQ  . (A-8) 

Figure  A-l  shows  plots  of  Eqs.  A-5  and  A-6  in  terms  of  SFL  and  arQ. 

Figure  A-l  may  be  used  in  the  following  way  to  assess  the 

general  importance  of  nonlinear  effects  on  the  propagation  of  spherical 

waves  in  air.  Nonlinear  effects  are  expected  to  be  weak,  moderate,  or 

strong  if  the  value  of  SFL  and  arQ  defines  a point  below,  between,  or 

above  the  two  curves,  respectively.  The  reason  is  this:  points  below 

the  two  curves  define  waves  for  which  r <r,  i.e.,  shocks  never  form; 

points  between  the  two  curves  depict  waves  for  which  r<r  <£,  i.e., 

max 

shocks  form  but  a sawtooth  does  not;  points  above  the  two  curves  define 

waves  for  which  r<r  , i.e.,  both  shock  and  sawtooth  formation  is 
max 

possible.  The  points  B and  C in  Fig.  A-l  represent  the  operating  points 
for  the  experiments  presented  in  Appendices  B and  C,  respectively,  and 
are  discussed  there. 

We  emphasize  here  that  Fig.  A-l  contains  no  Information 
regarding  the  length  scale  for  the  development  of  the  nonlinear  effects. 
That  is,  one  may,  for  example,  find  that  for  a given  set  of  source 
parameters  the  corresponding  values  of  SFL  and  ctro  define  a point  well 


FIGURE  A- 1 

DESIGN  CURVE  TO  ASSESS  THE 
IMPORTANCE  OF  NONLINEARITY  ON 
SPHERICAL  WAVE  PROPAGATION 

(CONDITIONS  ARE  FOR  AIR  AT  20°C) 
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above  the  r=r  curve.  Calculation  of  the  distances  r and  r may, 
max 

however,  show  that  r falls  well  beyond  the  distance  at  which  experimental 
observations  will  be  made.  Hence,  nonlinear  effects  would  be  of  little 
importance  for  the  given  length  scale.  In  summary,  if  use  of  Fig.  A-l 
shows  that  nonlinear  effects  may  be  important,  then  computation  of  the 
range  parameters  themselves  must  be  made  to  set  the  length  scale  of  the 
development  of  the  effects. 

In  most  practical  cases  the  acoustic  source  is  not  a vibrating 
sphere  but  a directional  radiator,  such  as  a baffled  piston  or  an  acoustic 
horn.  Though  we  have  considered  an  omnidirectional  source,  our  results 
may  be  extended  to  directional  waves  by  replacing  p1Q  in  Eq.  A-l  with 
P^qD(0),  where  D(0)  is  the  normalized  directivity  function.^0  But  the 
value  of  rQ,  the  effective  source  radius,  needs  to  be  specified.  If 
there  is  no  distortion  in  the  nearfield  of  the  radiator,  it  is  appropriate 
to  take  rQ= Rq,  where  is  the  length  of  the  nearfield.  Unfortunately, 
however,  this  is  seldom  the  case.  Nonlinear  effects  cause  the  signal  to 
distort  as  it  propagates  through  the  nearfield.  Several  researchers 
have  investigated  this  problem  empirically  and  found  that  choosing  rQ  in 
the  range  Fo/3<ro<3Ro/A  works  fairly  well  for  piston-type  radiators.^ 
Apparently  choosing  rQ  in  this  range  gives  about  the  right  weight  to  the 
distortion  that  takes  place  in  the  nearfield. 

At  this  point,  we  include  experimental  data  to  illustrate  the 
validity  of  using  the  three  critical  distances  to  assess  the  importance 
of  nonlinear  effects.  Unfortunately,  the  data  we  have  is  for  plane  waves, 
not  spherical  waves.  However,  the  basic  concept  is  the  same,  no  matter 
what  the  geometry  of  the  wave  motion  may  be.  Figure  A-2  shows  waveforms 
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MEASURED  WAVEFORMS  CORRESPONDING  TO  WEAK, 
MODERATE,  AND  STRONG  NONLINEARITY 
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for  the  case  of  weak,  moderate,  and  strong  nonlinearity.  The  condition 
Xmax^*  holds  for  the  first  set  of  wavesforms,  Fig.  A-2(a).  Under  this 
condition,  shocks  are  not  expected.  Notice,  indeed,  that  i shock  never 
forms.  The  steepening  at  x=x  is  less  than  that  measured  at  the  H°cond 
distance.  The  reason  is  that  the  high  frequency  components  damp  out 
before  a shock  can  form.  Figure  A-2(b)  shows  measured  waveforms  consist- 
ent with  the  condition  x<x  <x.  At  the  second  distance  (x=x)  considerable 
steepening  is  present  though  the  shock  is  not  fully  formed.  At  the  third 
distance  (x=x)  the  wave  is  slightly  more  distorted  than  at  the  second 
distance,  but  certainly  it  is  not  a full  fledged  sawtooth.  The  peculiar 
wave  shape  at  x=x  is  caused  by  dispersion  and  is  discussed  in  detail  in 
Chapter  2.  Figure  A-2(c)  depicts  the  case  x<x  . At  the  distance  x=x 
shock  formation  is  very  apparent,  and  the  characteristic  sawtooth  shape 
is  evident  at  x=x.  We  conclude  that  the  critical  distance  method  is  in 
good  qualitative  agreement  with  the  data  of  Fig.  A-2.* 

In  summary  Fig.  A-l,  used  in  conjunction  with  Eqs . A-2  through 
A-^,  should  be  quite  useful  in  estimating  the  importance  of  nonlinearity 
on  the  propagation  of  spherically  spreading  sound  waves.  While  the 
approach  is  not  new  (see  Ref.  3),  the  particular  graphical  representation 
is . 


Experimental  data  for  spherical  waves  is  given  in  Ref.  28. 
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APPENDIX  B 

A PERTURBATION  SOLUTION  OF  BURGERS'  EQUATION 
FOR  SPHERICAL  WAVES 


Our  purpose  here  is  to  obtain  a theoretical  prediction  valid 

for  weak  spherical  waves  in  a thermoviscous  fluid.  Several  other 

researchers  have  obtained  results  valid  for  the  weak-wave  problem. 51-33 
9 

Pernet  and  Payne,  for  example,  extended  the  phenomenological  approach  of 

23  xx 

Thuras  et  al.  to  the  spherical  wave  problem.  Blackstock  and  Willette, 

however,  obtained  a perturbation  solution  of  Burgers ' equation  valid  to 
third  order  in  the  source  Mach  number  e.  The  work  we  present  here  is  an 
extension  of  Blackstock' s et  al.  perturbation  analysis  to  fifth  order 
in  e.  We  also  show  a limited  test  of  the  theoretical  results  by  com- 
parison with  data  from  an  outdoor  propagation  experiment. 

The  Burgers'  equation  for  spherical  waves  is  (see,  for  example. 

Ref.  1) 

1 d(ru)  a d2u  p du  /0  ,v 

7*37“  - -57-2  = • (B-1) 

to  at  c 

o 

where  a is  the  thermoviscous  absorption  coefficient  at  angular  frequency  05 
and  t'=t-(r-r 0/cQ)  is  the  retarded  time.  Let  the  boundary  condition  be 

u(rQ,t)  = u1Q  sin  tofc  . (B-2) 
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It  is  convenient  here  to  cast  Eqs.  B-l  and  B-2  in  terms  of  the 
dimensionless  quantities  V=(r/rQ) (u/cq)  and  y=ojt'.  Jfaking  the  substitu- 
tions for  V and  y and  also  defining  5=pkrQ,  Eqs.  B-l  and  B-2  become 


dV 

3? 


a 


—4-  = ivr 

r oy 


dv 


(B-3) 


and 


V(rQ,y)  = € sin  y , (B-U) 

respectively.  We  seek  a solution  of  Eq.  B-3  in  the  form  of  a perturbation 
series.  That  is,  we  suppose  a solution  of  the  form 


V . i e\<">  . (B-5) 

n=l 

Substituting  Eq.  B-5  into  Eq.  B-3,  we  obtain 


00 


E 

n=l 


a 


sV"> 

T7~ 


.i 

r 


E 

i+j=n 


r(  1 ) 


(d  )' 

y 


= o 


(B-6) 


In  order  that  Eq.  B-6  be  satisfied,  the  bracketed  term  must  vanish  for  each 
value  of  n.  In  general  the  resulting  differential  equations  corresponding 
to  terms  of  order  en  (denoted  0(en))  may  be  written  as 


0(€D) 


dv(n) 

+ “3T" 


- 1 IT 

6y2  r i+j=n 


v(i)  Sv(j) 
dy 


(B-7) 


We  see  that  in  general  the  equations  must  be  solved  in  sequence.  For 

(2) 

example,  the  source  terms  for  the  second  order  solution  V'  depend  on 
V^.  The  source  terms  for  the  third  order  solution  depend  on 

and  . We  now  proceed  to  solve  the  first  five  such  equations. 
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The  first  order  equation  is 


d^y(l) 


- a 


A) 


= 0 


(B-8) 


Hence,  the  first  order  solution  Vv  ' , which  satisfies  the  boundary 
condition  Eq.  B-U,  is 

/ v -a(r-r  ) 

V = e ° sin  y , (B-8a) 

which  is  the  familiar  solution  of  linear  acoustics.  The  second  order 


equation  is 


g - L y(D  & 


(1) 


-L- 

2r 


*T 

-2a(r-r  ) 


(B-9) 


sin  2y 


,(2). 


The  solution  of  Eq.  B-9  may  be  found  by  assuming  Vv  =F(r)  sin  2y. 
The  result  is 


v(2)  . -Mr-p0)  t 


22  2y  , 


(B-9a) 


where 


2a( r ' -r  ) 
o 


22 


l:--r 

» n 


dr' 


We  shall  have  occasion  to  write  down  several  integrals  of  the 
form  I and  adopt  the  following  convention.  The  Integral  I means 

ran  mn 

the  Integral  which  must  be  evaluated  to  obtain  the  contribution  to  the 
nth  harmonic  from  the  mth  order  perturbation. 
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. a d2v(U)  = £ [>)  ay(g)  a(v(1)v(3))'[ 

dr  sy2  r |_  ~3y~  3y  J 


% [(’ 


= £ Il'V6 


-2a( r-r  ) 
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• 5I33  6 


-10a(r-r  ) 


) 


sin  2y 


(' 


Mx22  e 


-8a(r-r  ) 


+ 6l33  e 


-10a( r-r  ) 


) 


sin  4y 


(B-ll) 


,(M  . £ [. 


-Mr-r  ) 

(i42+3i'i+2)  e 0 sin  2y 


-l6a(r-r  ) 


+ (144+  6l  \k)  e u sin  4y 

2a(r'-r  ) 


] 


k2 


I. 


r I51(r’)  e 


dr' 


(B-lla) 


-6a(r'-r  ) 

rl  (r’)e 

dr' 


L44 


■/. 

■/,: 


P 8a(r'-r  ) 

2^i\  „ o' 


r I22  (r')  e 


dr' 


:\u  “ ft 


6a(r'-r  ) 

r I,,(r' ) e 


dr’ 


Because  of  the  rather  formidable  algebra  involved,  the  only 

part  of  the  fifth  order  solution  that  was  calculated  was  the  contribution 

( 5) 

to  the  fundamental.  This  contribution  is  denoted  Vi.  ' . . The  results 

fund. 

are  listed  below. 
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a ^ i a(v^V4))  + a(y(2V3)) 


2 r 


, ■ s'  •*'’•’  [■„  - J ■•»] 

[ T I22^ 

/ I42(r’)  + 5I'42(r’)  + 

^ ro  L 


>•>  yr'>1 

2 I r’ 


r Igg (r1)  I^(r' ) e 


-12a(r'-r  ) 


As  may  toe  seen  from  the  higher  order  results,  the  solutions  are  easy  to 

write  down  in  terms  of  the  integrals  T . It  is  the  calculation  of  these 

mn 

integrals  which  is  indeed  a formidable  task. 

Let  us  now  put  our  results  in  terms  of  the  normalized  harmonic 
amplitudes  Bn= ( r/ r0 ) ( un/ u0 ) * Defining  ao=e£,  we  obtain 


-a(r-r  ) r a ^ a ^ 

Bi  = e ° I1  ” TP  x3l  + “TT 


§"  (T51  ' 2 I'5l)J 


-^(r-r  )pJ  a * | 

B2  " * ° h2  J22  " T (I42*5I'‘t2J 


B5  = 6 


-90(r-r0)  r5ao 

T-  I5jJ  ’ 
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-l6ot(r-r  ) [a  3 1 

B4  ~ e ° Ijx-  ' (B-16) 

As  a check  on  our  theoretical  results,  we  now  inquire  as  to 

the  behavior  of  the  solutions  B in  the  limit  of  small  values  of  ar. 

n 

We  expect  that  these  limiting  forms  should  correspond  to  the  harmonic 
amplitudes  found  from  the  Bessel-Fubini  solution  (see,  for  example, 

Ref.  12),  which  are  given  by 


B„  - h Jn  <■">  - <B‘1T) 

where  a=0ekroln(r/ro) . To  compare  our  solution  to  Eq.  B-17,  we  first 

establish  the  limiting  forms  of  the  integrals  I . In  the  limit  as 

mn 

ar-*0,  the  integrals  I reduce  to  the  following. 


I22  = ln  r 
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t * 2 r 

• ln  r; 


*33  ■ 5 l”2  'o 


T 1 ln3  r 

zk2 = z 1x1  r 
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1 ,.3  £. 
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Vk2  = 


T . I L 
r44  3 In  r 


I'w.-51"5r 


(B-l8a) 
(B-l8b) 
(B-l8c) 
(B-l8d) 
(B-l8e) 
(B-l 8f) 


(B-l8g) 
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J51  = S ^ r (B-l8h) 

o 

I,5i  = HlnUr  • (e-ieo 

o 

The  expressions  for  the  harmonic  amplitudes  B are.  in  the  limit  as  ar-»0  , 

n ' 


2 k 

Bi  = 1 • T + W ’ (B-19) 

3 

B2  = 2 ' T ’ (B-20) 


(B-21) 
( B-22) 

Comparison  of  Eqs.  B-19  through  B-22  with  the  series  expansion  of  Eq.  B-17 
shows  that  our  results  are  consistent  with  the  Bessel-Fubini  solution  to 
the  orders  of  a which  we  have  kept.  For  example,  Eq.  B-19  is  the  series 
expansion  of  Eq.  B-17  to  fourth  order  in  a.  Hence,  the  results  of  the 
perturbation  solution  (Eqs.  B-1J  through  B-l6)  do  indeed  reduce  to  the 
Bessel-Fubini  solution  as  Ctr-»0 . 

In  many  cases  where  such  a perturbation  solution  is  useful, 
some  distortion  is  evident  in  the  source  waveform  V(rQ,y).  Commonly 
this  distortion  is  primarily  second  harmonic  (especially  if  it  is  due  to 
nonlinear  propagation  distortion  in  the  nearfield).  We  shall  now 


B - 50 
B3  ~ “8" 


*-r 
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determine  how  the  presence  of  an  "initial"  second  harmonic  signal  changes 
the  lower  order  perturbation  results.  Suppose  that,  in  place  of  Eq.  B-4, 
the  boundary  condition  is 


V(rQ,y)  = e 


sin  y + a sin 


(2y+<p)j 


(B-23) 


where  a and  cp  are  the  relative  amplitude  and  phase  of  the  second  harmonic, 
respectively.  The  first  order  perturbation  solution  satisfying  this 
boundary  condition  is 

m <t(r-r  ) -4a(r-r  ) 

V'  ' = e sin  y + a e sin(2y+cp)  . (B-24) 

(2) 

The  second  order  solution  V has  a source  term  proportional  to 
V^^dV^V^y).  Hence,  it  is  easily  seen  that  V ^ contains  contributions 
to  the  first,  second,  third,  and  fourth  harmonics.  The  expression  for 


at  ■0,(r*po)  t -Mr-r  ) 

* 2 * *21  sin(y+!p)  + * e Ig2  sin  2y 


+ P(r)  sin(3y+q>)  + G(r)  sin(4y+2qp) 


(B-25) 


where 


*21 


-4a(r'-ro) 


dr' 


(2) 


We  note  here  that  the  second  harmonic  component  of  V 
the  "initial"  second  harmonic.  The  expressions  for  B^ 


is  not  affected  by 
and  Bg  are,  then, 


to  second  order 
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and 


(B-26) 


(B-27) 


We  now  vlsb  to  compare  the  results  of  the  perturbation  analysis 

with  data  from  an  outdoor  propagation  experiment.  In  particular  we  wish 

to  compare  the  measured  values  of  p^  and  Pg  with  the  perturbation  analysis 

for  a boundary  condition  similar  to  Eq.  B-23.  The  source  was  an  array  of 

p A 

seven  electroacoustic  drivers  with  exponential  horns.  A side  view  of 
the  array  Is  shown  In  Tig.  B-l.  The  diameter  and  transmitting  frequency 
of  the  array  were  0.5  m and  8.25  kHz,  respectively.  The  absorption 
coefficient  a was  O.OO76  Np/m.  The  amplitudes  p^  and  pg  at  the  beginning 
of  the  far fie Id  (rQ=6.1  m)  were  551  nbar  and  1^3  pbar,  respectively, 
and  the  phase  angle  ip  was  observed  to  be  zero.  In  terms  of  the  parameters 
discussed  In  Appendix  A,  the  operating  conditions  for  this  experiment  are 
SFIa=159»8  dB,  otro=0.046.  These  conditions  determine  a point  on  Pig.  A-l 
(point  B)  that  falls  below  the  Ts‘rmx  curve;  therefore  the  perturbation 
analysis  should  be  valid.  Figure  B-2  shows  the  measured  data  and  theoreti- 
cal results.  Since  the  wave  Is  weak,  the  linear  theory  result  (Eq.  B-8e) 
describes  the  behavior  of  p^  fairly  well,  though  the  second  order 
perturbation  result  (dashed  curve)  provides  a slightly  better  fit  at  the 
greater  ranges.  The  lower  solid  curve  Is  the  linear  theory  prediction 
for  Pg.  The  dotted  curve  Is  the  second  order  perturbation  result  for  pg 
for  no  Initial  second  harmonic  (Eq.  B-9a).  Since  the  phase  dl  Terence  cp 
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FIGURE  B-l 

ACOUSTIC  ARRAY  USED  FOR 
WEAK-WAVE  PROPAGATION 


Taken  from  Ref.  28 
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Is  zero,  the  correct  second  order  perturbation  result  (dashed  curve)  for 
Pg  is  the  arithmetic  sum  of  the  dotted  curve  and  the  lower  solid  curve. 
The  need  for  taking  account  of  the  initial  second  harmonic  distortion  is 
well  illustrated.  It  is  also  interesting  that  although  the  fundamental 
is  well  described  by  linear  theory,  the  linear  theory  prediction  for  the 
second  harmonic  is  10  d£  below  the  data  at  the  greatest  range.  The 
overall  agreement  as  regards  the  predicted  and  measured  second  harmonic 
is  excellent. 

In  conclusion  we  have  presented  a fifth  order  perturbation 
solution  of  Burgers'  equation  for  spherically  spreading  sound  waves. 
Though  the  integrals  involved  become  increasingly  more  complicated  for 
the  higher  order  results,  the  lower  order  results  should  be  useful  for  a 
great  many  weak-wave  problems. 


APPENDIX  C 

SOLUTION  OF  THE  RUDNICK  EQUATION  FOR  SPHERICAL  WAVES 


34 

Blackstock  suggested  a spherical  wave  version  of  the  Rudnick 
decay  rate  equation  (Eq.  3-7).  By  analogy  with  the  plane  wave  case, 
Blackstock  assumed  that  the  decay  rate  dp^/dr  of  the  fundamental  pressure 
component  is  the  sum  of  the  decay  rate  due  to  absorption  (dp ,/dr)  = -cep 

X cLDS  • X 

and  the  decay  rate  (dp1/dr)„  of  a spherically  spreading  sound  wave  of 
finite  amplitude.  The  latter  decay  rate  may  be  found  as  follows.  Let 
the  boundary  condition  be 

P1(rQ,t)  = p1Q  sin  cut  . (C-l) 

The  value  of  in  the  sawtooth  region  (0*3)  in  the  absence  of  absorption 
is 

(C-2) 
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The  decay  rate  (dp./dr).  is,  from  Eq.  C-2 
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8o 


By  supposition  the  total  decay  rate  dp^/dr  is  then 


f .a. 


_ 1 pek  2 

- ? pi ' pi 


(c-M 


We  now  wish  to  solve  Eq.  C-1*  subject  to  the  condition  that  the  solution 
reduces  to  the  weak  shock  solution  (Eq.  C-2)  as  a-»0  . The  form  of  the 
solution  for  the  plane  wave  case  (Eq.  3-11 ) suggests  that  we  try  a solution 
of  the  form 


Pt  = — 


-a(r-r  ) 
r . o 

o A e 


T V f(r) 


(C-5) 


where  A is  a constant  and  the  function  f(r)  is  to  be  determined. 
Substitution  of  Eq.  C-5  into  Eq.  C-4  leads  to  the  following  form  of  f(r): 


f(r) 


A pekr 
^10 


(C-6) 


o 

As  in  the  plane  wave  case,  the  choice  of  A=2p^Q  yields  the  weak-shock 
solution,  Eq.  C-2,  as  qmO.  The  definition  of  the  exponential  integral 
E^a)  (see,  for  example,  Ref.  35), 


dX 


may  be  used  to  write  the  expression  for  f(r)  as  fo^ows: 


Or 

f(r)  = pekro  e ° [E^Or^-E^otr)  ] . (C-7) 
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The  expression  for  p^,  found  from  Eqs.  C-5  and  C-7>  is 


r 

o 

Pl  * T 


^10  e 
otr 


-a(r-r  ) 
o 


(C-8) 


1 + P€krQ  e °[E1(aro)-E1(ar)  ] 

Equation  C-8  is  simple  to  evaluate  in  practice.  The  exponential 
integral  E^Ofr)  may  be  read  from  tables  or  evaluated  by  certain  series 
representations.  See,  for  example,  Ref.  35.  Indeed,  in  the  range 
O^arSl  the  series  representation  for  E^Cfr)  may  be  written  as^ 


E (Or)  = -ln(otr)  + £ a (ocr)1 
i=0 


(C-9) 


with  a maximum  percentage  error  of  9x10  . 

The  proposed  solution  (Eq.  C-8)  was  compared  to  freefield 

propagation  data  taken  using  a siren  operating  at  6.1  kHz  in  air.  Two  views 

of  the  siren  are  shown  in  Fig.  C-l.  The  20  circular  ports  in  view  b are 

the  sound  sources.  The  siren  operates  by  periodically  interrupting  the 

air  flow  through  the  ports  by  means  of  a motor-driven  titanium  steel 

rotor  located  just  behind  the  ports.  For  further  details  see  Ref.  28. 

Because  of  the  relatively  high  acoustic  power  output  of  the  siren 

(approximately  600  W of  acoustic  power  with  a 40°  beamwidth  between 

the  3 dB  down  points),  the  characteristic  sawtooth  wave  shape  was  already 

in  evidence  at  a propagation  distance  of  approximately  2 m.  Our  measured 

boundary  condition  was,  therefore,  a sawtooth  wave  at  a given  distance  r^ 

from  the  source.  To  use  Eq.  C-8  we  must  compute  from  our  measured  data 

the  amplitude  p1Q  and  effective  source  radius  rQ  of  a sinusoid  which 

would  yield  the  measured  sawtooth  wave  at  range  r . This  is  easily 

m 
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FIGURE  C-l 

ACOUSTIC  SIREN  USED  IN  THE 
STRONG-WAVE  EXPERIMENT 
(o:  SIDE  VIEW,  b:  VIEW  FROM  ABOVE  THE  SIREN) 


Token  from  Rof.  28 
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accomplished  by  using  the  sawtooth  solution  (Eq.  C-2).  Let  p,  be  the 

lm 

measured  value  of  the  fundamental  pressure  component  at  range  r . Then, 

m y 

using  Eq.  C-2,  we  have 


and 


a = 3 

r=r 

m 


a 


Pp 


lO^o 


P c 
*o  o 


(C-10) 


(C-ll) 


Equations  C-10  and  C-ll  may  be  solved  for  p^q  and  rQ,  since  all  other 
quantities  are  known. 

Figure  C-2  shows  the  results  of  a comparison  of  Eq.  C-8  with  the 

measured  data.  The  value  of  the  absorption  coefficient  a was  found  to  be 

0.0059  Np/m  from  a small-signal  experiment  under  similar  conditions  and  is 

also  consistent  with  the  computed  value  from  Ref.  36.  The  values  of  p, 

lm 

and  were  2.2^C10^  pbar  and  2.19  m>  respectively.  Hie  values  of  p1Q  and 
rQ  computed  from  Eq.  C-10  were,  respectively,  1.1x10^  pbar  and  0.09  m. 

The  solid  curve  is  the  linear  theory  prediction;  the  dashed  curve  represents 
the  solution  of  the  "Rudnick  model"  equation.  The  latter  is  seen  to 
provide  a good  fit  to  the  data. 

In  summary  the  model  proposed  here,  although  ad  hoc,  does  seem 
to  provide  a good  description  of  the  experimental  data.  Further  experi- 
mental verification  of  this  model  may  be  found  in  Ref.  28.  The  model 
should  be  useful  in  a number  of  spherical  wave  problems,  though  its  use 
is  restricted  to  strong  waves. 
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FIGURE  C-2 

COMPARISON  OF  THE  RUDNICK  MODEL 
WITH  PROPAGATION  DATA 
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